UAH Technical Report No: 

5-33020 

Award No: 

NAS8-38609 

Delivery Order No: 

49 


DESIGN OF A NEW HIGH-PERFORMANCE 
POINTING CONTROLLER 
FOR THE HUBBLE SPACE TELESCOPE 


Final Report 

3 August, 1992 through 3 October, 1993 


Prepared by 
Dr. C. D. Johnson 

Electrical and Computer Engineering Dept 
University of Alabama in Huntsville 
Huntsville , Alabama 35889 


Prepared for 

Precision Pointing Systems Branch 
NASA, Geo. C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


October 1993 



NASA 

Niux w ni i io yc 

^Dct«uniw*tr 


Report Documentation Page 


1. Report No. 


2. Government Accession No. 


4. Title and Subtitle 


Design of a New High-Performance Pointing Controller 
for the Hubble Space Telescope 


3. Recipient s Catalog No. 


5. Report Date 

October 1993 


6. Performing Organization Coae 


7. Authorial 


Dr. C. D. Johnson 


9. Performing Organization Name and Address 

University of Alabama in Huntsville, ECE Dept. 
Huntsville, Alabama 35899 


12- .Sponsoring Agency Name and Address 


National Aeronautics and Space Administration 

Washington, DC 20546-0001 

Marshall Space Flight Center, Huntsville, AL 35812 


8. Performing Organization fteoort No. 


10. Work Unit No. 

D.O. #49 


11. Contract or Gram No. 


NAS8-38609 


13. Type of Report and Period Covered 

FINAL: 8/3/92-10/3/93 


14. Sponsoring Agency Code 


IS. S up pl emen tary Notes 


16. Abstract 


A new form of high-performance, disturbance-adaptive pointing controller for the 
Hubble Space Telescope (HST) is proposed. This new controller is all linear (constant 
gains) and can maintain accurate "pointing" of the HST telescope in the face of persistent, 
randomly triggered uncertain, unmeasurable "flapping" motions of the large attached solar 
array panels. Similar disturbances associated with antennas and other flexible appendages 
can also be accommodated. The effectiveness and practicality of the proposed new 
controller is demonstrated by a detailed design and simulation testing of one such 
controller for a planar-motion, fully nonlinear model of HST. The simulation results show 
a high degree of disturbance isolation and pointing stability. 


17. Key Words (Suggested by Autnoris)) j 18. Distribution Statement 

Hubble Space Telescope j 

Pointing Controller for HST i 

Pointing Control System 

j 

j 

1 19. Secunty Gassif. (of inis report) I 20. Security Classif. lot tms page) j 21. No. of pages ; 22. Price 


NASA FORM 1626 OCT 86 








DISCLAIMER 


The views and opinions expressed in this report are those of the 
Principal Investigator and his collaborator, Mr. Stewart Addington, and are 
not necessarily the views and opinions of NASA, in general, or the Precision 
Pointing Systems Branch of Marshall Space Flight Center, in particular. 



OVERVIEW OF THIS REPORT 


This report describes a new methodology for designing a disturbance- 
adaptive, precision pointing controller for the Hubble Space Telescope 
(HST). The proposed new controller can maintain a consistent high-q ualit y 
of pointing stability in the face of persistent, randomly triggered, "flapping" of 
the two large flexible solar arrays attached to either side of HST. Similar 
HST unknown, unmeasurable "disturbances" associated with deflections of 
other flex-body appendages (long, slender antenna booms, etc.) can also be 
accommodated by the same pointing controller. 

A noteworthy feature of this new disturbance-adaptive pointing 
controller is that the structure of the controller algorithm is all linear, has all 
constant coefficients and is relatively low-order, thus enhancing the 
controller's reliability and implementation attributes. 

The effectiveness of the proposed new precision pointing control 
methodology is demonstrated in this report by a detailed design, and closed- 
loop simulation testing, of one such pointing controller for a planar-motion 
(single-axis), "exact" non-linear model of HST. The simulation results show 
a high-degree of pointing stability in the face of persistent flex-body 
"flapping" motions of the simulated solar arrays. A novel feature of the 
proposed controller is its dual-mode capability. In particular, the controller 
has the capability of switching into a special non-pointing control mode in 
which the pointing control actions "sway" the telescope’s main-body in a 
series of back-and-forth motions that are automatically timed and orchestrated 
to induce "active damping" (damping augmentation) of the flapping solar 
arrays. This action results in a significant and rapid reduction of the solar 
array oscillations. 



The rationale, general design procedure, and worked (single-axis) 
example for the proposed new dual-mode precision pointing controller is 
presented in Chapter 1 of this Report. Details of the derivation of the "exact,” 
non-linear HST model (planar-motion) and the associated model equations 
are presented in the various Appendices of Chapter 1, along with various 
design formulae. 

The scientific basis of the precision pointing controller design 
methodology proposed herein is relatively new and/or unfamiliar to some 
practicing control designers. Therefore, to aid such readers in becoming more 
familiar with the underlying scientific basis of our methodology, and to 
explain more fully why this new approach is essentially the "best" one can do 
in designing precision pointing/tracking controllers for systems operating in 
uncertain disturbance environments, a tutorial account of the control-theoretic 
principles of our methodology is presented in Chapter 2. 

The detailed design and simulation results presented in this report were 
developed in collaboration with Mr. Stewart Addington, Graduate Student in 
the ECE Dept, of UAH; see also the Acknowledgment at the end of this 
chapter. 
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Chanter 1 

DESIGN OF A NEW HIGH-PERFORMANCE 
POINTING CONTROLLER 
FOR THE HUBBLE SPACE TELESCOPE* 

Chapter Summary 

Cyclic thermal expansions and mechanical stiction effects in the 
Solar Arrays on the Hubble Space Telescope (HST) are triggering repeated 
occurrences of damped, relaxation-type flex-body vibrations of the solar 
arrays. Those solar array vibrations are, in turn, causing unwanted, 
oscillating disturbance torques on the HST main body, which cause unwanted 
deviations of the telescope from its specified pointing direction. 

In this chapter we propose two strategies one can adopt in 
designing a telescope-pointing controller to cope with the aforementioned 
disturbances: (i) a "total isolation" (TI) control strategy whereby the HST 
controller torques are designed to adaptively counteract and cancel-out the 
persistent disturbing torques that are causing the unwanted telescope motions, 
and (ii) an "array damping" (AD) control strategy whereby the HST 
controller torques are used to actively augment the natural dampening of the 
solar array vibrations and the attendant telescope motions, between 
triggerings of the stiction-related flex-body relaxation oscillations. 

Using the principles of Disturbance-Accommodating Control 
(DAC) Theory a dual-mode pointing controller for a generic, planar-motion 

* The control design and simulation results presented in this chapter were developed in collaboration with 
Mr. Stewart Addington, Graduate student in the ECE Dept of UAH; see also the Acknowledgment at the 
end of this chapter. 


UAH / ECE Dept / Dr. C. D. Johnson 


1-1 


Chapter 1, Final Report 




Design of a New High-Performance Pointing Controller for the Hubble Space Telescope 

NAS8-38609/D. 0.49 


(single-axis) model of the HST is proposed. This pointing controller 
incorporates both the TI and AD modes of disturbance-accommodation. 
Simulation studies of the closed-loop system using generic parameter values 
clearly indicate, qualitatively, the enhanced pointing-performance such a 
controller can achieve. 

1. INTRODUCTION 

The two large deployable solar array panels attached, in a cantilever- 
beam fashion, to opposite sides of the Hubble Space Telescope (HST) main 
body (Figure 1) have introduced a new dimension to the problem of 
controlling the HST's pointing direction. In particular, the uneven thermal 
expansions of the solar array's collapsible structural members, caused by 
cyclic solar heating effects, and mechanical coulomb friction/stiction effects 
that resist thermal expansions between structural members are causing 
persistent, jerky, "flapping motions" of those solar arrays. Those flex-body 
flapping motions are of a damped-oscillation type and are triggered in a 
sporadic, random-like manner, depending on the intensity of the solar 
heating/cooling and on the varying stiction thresholds associated with sliding 
interactions between various structural members of the arrays. 

These flapping motions of the solar arrays induce, through their 
attachment points with the main body of the HST, a series of randomly 
triggered damped-oscillation disturbance torques that cause the HST main 
body to veer away from the precision point needed to satisfy experiment 
requirements. Thus, an effective HST pointing control system must cope with 
this uncertain disturbance environment. 

In this chapter we will derive a new conceptual form of HST pointing 
controller that effectively accomplishes this goal. Simulation results, using a 
non-linear, planar-motion model of the HST with generic parameter values. 
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Fig. 1 : The Hubble Space Telescope (HST) with Deployed Solar Arrays 
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demonstrate the effectiveness of the proposed new controller in a single-axis 
control context. A generalization of the control concepts presented here, to 
include full 3-axis pointing control of the HST, should provide a significant 
enhancement to the actual HST pointing performance. 

2. CHARACTERIZATION OF THE SOLAR ARRAY FLAPPING 
DISTURBANCES IN THE HST CONTROL PROBLEM 

The randomly-triggered flapping solar arrays pose a unique challenge 
to the HST control system design because the uncertain disturbance torques 
they induce to the HST main body are not effectively modeled as steady-state 
"random processes" with known means and variances. Consequently, 
traditional "stochastic control theories" [1] are ineffective in accomplishing 
high-performance, precision pointing of HST in the face of those flapping 
solar arrays. 

The essential features of the HST uncertain disturbance torques due to 
flapping solar arrays are: (i) once triggered they have a distinguishable 
"damped-oscillation" type waveform characteristic that may have varying 
initial conditions, but is always relatively smooth and well-behaved 
(compared to the grossly jagged, erratic "white noise" of random process 
theory), and (ii) they are triggered in a once-in-a-while fashion with a 
relatively wide time-window between successive triggers (compared to the 
virtually continuous fusillade of triggering events associated with 
mathematical "white noise"). Thus, for instance, individual components 
(modes) wj(t) of the HST uncertain disturbance torques, due to the solar 
arrays, can be rather accurately represented (modeled) by a mathematical 

"spline" expression of the form 

Wj(t) = Ci e - ^ t sin(co i t) + C 2 e - ^ t cos(Oj t) (1) 
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where (3j>0 denotes the damping parameter and ©jX) denotes the oscillation 
frequency and where the spline "weighting-constants", C 1? C 2 in Eq.(l) are 

allowed to jump in value from time-to-time, in an uncertain, random-like 
manner (thereby allowing for different "initial conditions" on Eq.(l) at each 
triggering event). The numerical values of (pi,©0 in Eq.(l) can be estimated 
from examination of actual flight data, results of ground experiments, etc. 

It can be argued that, in practice, the parameters pj,©i in Eq.(l) might 
be subject to some degree of uncertain, time-dependent variation due to wear, 
distortions, thermal effects, etc. which may render the model Eq.(l) 
ineffective. In such cases, an effective alternative to Eq.(l) is the 
considerably more robust polynomial-spline representation 

Wi(t) = C} + C 2 t + C 3 t2 + ... + Cm 1 (m * 1 > (2) 

where m is an appropriate positive integer and, as in Eq.(l), the "weighting- 
constants" Ci, C 2 ,...,C m are allowed to jump in value in an uncertain, once-in- 
a-while fashion, hereafter called "stepwise-constant" behavior. Expression 
Eq.(2) can effectively represent a broad class of uncertain, meandering-type 
disturbance functions wj(t), including complex "oscillations" of the type 
Eq.(l) with completely unknown waveform characteristics. In practice, 
satisfactory results using Eq.(2) can usually be achieved even with relatively 
small values for m, typically m = 3 or m = 4 (the so-called quadratic and 
cubic "spline models", respectively.) 

In summary, individual modes of the uncertain disturbance torques 
associated with the flapping solar arrays on the HST have relatively smooth, 
well-behaved, time-domain waveform characteristics of the randomly 
triggered, damped-oscillation type. Consequently, those disturbance torques 
are effectively represented by mathematical "spline expressions" of the form 
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Eq.(l) or Eq.(2). It is remarked that in some situations, a combination of 
Eq.(l) and Eq.(2) may be most effective. 

3. ACCOMMODATION OF SOLAR ARRAY DISTURBANCES IN THE 
HST 

As stated earlier, the HST disturbance torques associated with flapping 
solar arrays cause the point of the HST main body to deviate from the desired 
direction. Thus, it is clear that those persistent disturbances cause only 
unwanted, upsetting effects. Consequently, the primary goal of the HST 
pointing controller, with respect to accommodating those disturbances, should 
be to (ideally) generate opposing control torques that automatically adapt-to 
and cancel-out (counteract) the persistent disturbance torques and their 
upsetting effects, in real-time. Under this latter mode of control, the main 
body of the HST would be effectively isolated from the disturbing motions of 
the flapping solar arrays, so that the point of the HST could then be regulated 
as if there were no solar array flapping motions. Hence, we hereafter refer to 
this control strategy as the total isolation (TI) mode of control. 

In the TI mode of control, no effort is made to mitigate the flapping 
motions of the solar arrays themselves. However, the (mild) natural structural 
damping effects in the solar array structures will tend to dampen-out those 
oscillations— until such time that the oscillations are once again triggered by 
another thermal/stiction relaxation effect. In some situations, it may be 
desirable to employ the HST controller to hasten the natural damping-out of 
the solar array oscillations. This can be accomplished by designing the HST 
controller to create strategic rocking motions of the HST main body that are 
so timed and phased as to accomplish ’’active damping” (damping 
augmentation) of the solar array oscillations. We hereafter call this control 
strategy the array damping (AD) mode of control. Practical realization of the 
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AD mode of control would seem to require a relatively "smart" type of HST 
controller. However, it will be shown that a simple linear controller will do 
the job. 

The IT and AD modes of control constitute what we believe to be the 
two main options to be considered in designing the HST pointing controller to 
accommodate uncertain disturbances such as oscillations of the solar arrays. 
For maximum flexibility in implementations, it is desirable to have a mode- 
switching arrangement for the TI7AD control laws, whereby one can 
gracefully shift control action from the IT mode to the AD mode (and vice- 
versa) in real-time, as dictated by the real-time needs on-station. 

In the subsequent sections of this chapter we will demonstrate how a 
unique branch of modem control (Disturbance-Accommodation Control 
[DAC] Theory [2]-[6]) enables one to systematically design physically 
realizable, all linear, time-invariant controllers that embody the TI and AD 
modes of disturbance-accommodation described in this section. 

4. A SIMPLIFIED, PLANAR-MOTION CONFIGURATION MODEL 
OF THE HST 

The objective of this chapter is to develop, and demonstrate the 
performance characteristics of a new control concept for the HST pointing 
controller design. For this purpose, we will consider a simplified planar- 
motion (single-axis) generic configuration model of the HST, Figure 2, in 
which the inertias and flex-body (flapping) motions of the attached solar 
arrays are represented by rigid, movable arms attached to either side of the 
HST main body through idealized "pin-joints." To replicate one mode of the 
structural flexibility and structural damping effects of each of the solar arrays, 
the back-and-forth oscillations 0j(t) 02 (t) of the two arms about their 
respective pin-joints are considered to be resisted by linear torsional "spring" 
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Fig. 2: Planar-Motion (Single-Axis) Configuration Model of HST and Attached Solar Arrays 
(Replicating One Flex-Body Mode for Each Solar Array) 
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and linear viscous "damping" effects that may be different for each arm. The 
normal (equilibrium) positions of the arms are perpendicular to the HST main 
body. 

To simultaneously replicate multiple flex-bodv modes of the solar 
arrays, one simply imagines multiple rigid arms attached to each of the two 
pin-joints in Figure 2, where the rotational spring constant damping mass and 
inertia of each arm is chosen to conform to the known frequency, damping 
and inertial properties of the particular flex-body mode being replicated. In 
this way one can consider as many simultaneous flex-body solar array modes 
(and flex-body modes of other HST appendages) as desired, using the pin- 
jointed rigid-arm technique as shown in Figure 3. Our decision to consider 
here only one flex-body mode. Figure 2, for each solar-array is motivated by 
the concept-demonstration nature of this chapter. In particular, we wish to 
avoid introducing additional complexity that is not essential to illustrating the 
basic concept of our modeling and control-design. 

The planar-motion, concept-demonstration model in Figure 2 can be 
viewed as representing the actual HST when the solar arrays are rotated such 
that the flapping motions of each solar array consist of essentially one flex- 
body mode (which may be different for each array) and such that the 
associated disturbance torques induced to the HST main body cause relatively 
little out-of-plane rotational motions of the main body. The new controller 
scheme we develop here for the single-axis model in Figure 2 can be 
employed in each of the other two (rotational) axes (with proper coordination 
imposed between each such controller) to achieve full three-axis coordinated 
control of the HST main-body rotations for arbitrary positioning of the solar 
array panels. 
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Fig. 3: Planar-Motion (Single -Axis) Configuration Model of HST and Attached Solar Arrays 
(Replicating 3 Flex-Body Modes for Each Solar Array) 
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It is remarked that our introduction of the simplified, planar-motion 
HST configuration model in Figure 2 was not motivated primarily by a desire 
to keep things simple but, rather, by the technological necessity of having a 
mathematical model of HST that is based on "first principles" of dynamics. 
In particular, the ability to derive a controller that accomplishes the TI mode 
of disturbance-accommodation, as described in the previous section, requires 
that one begin with an HST mathematical model that embodies the actual 
dynamical equations of motion for the HST main body and for the "interface 
dynamics" associated with each of the attached solar arrays. Alternative 
mathematical models based on modal decompositions and/or transfer-function 
methodologies may not provide such detailed information and may "obscure" 
some novel control possibilities such as the TI mode of control we develop 
here. 

5. EXACT EQUATIONS OF MOTION FOR THE PLANAR-MOTION 
HST MODEL IN FIGURE 2 

The simplified, planar-motion HST configuration model in Figure 2 is, 
in fact, rather "complicated" from the dynamics point of view. In fact, the 
derivation of the exact equations of motion for Figure 2 would be extremely 
difficult and time-consuming using any of the classical methods of dynamics, 
such as Newton’s, Lagrange’s, or Hamilton’s methods. For this study, the 
exact equations of motion for Figure 2 were derived, in completely general, 
explicit symbolic form, using the method of Kane’s Equations [7] as 
implemented in the computer-aided modeling program Autolev©. The 
resulting mathematical model consists of a set of five (5) simultaneous, 
second-order, highly nonlinear, differential equations which can be written in 
the vector-matrix form 

M(y)y + Dy+ Ky + f(y,y) = bu (3) 
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where y is the 5-vector of position variables y = (<j>, 0 1? 0 2 , £, r|) as defined in 
Figure 2 and u is the (single-axis) HST main body control torque relevant to 
the planar-motion configuration shown in Figure 2. The elements of the 5x5 
matrices (M(y), D, K} and of the vectors (b,f(y,y)} in Eq.(3) are determined 

by various masses, inertias, lengths, etc., and their explicit expressions are 
given in Appendix A of this chapter. Those expressions are completely 
general in the sense that they are given in terms of arbitrary, symbolic masses, 
inertias, lengths, etc. 

The natural dynamic behavior of the mathematical model Eq.(3) turns 
out to be surprisingly complicated, rich in diversity, and often counter- 
intuitive, owing to the strong, non-linear "inertia-coupling” that exists 
between the motions {(J>(t), 0j(t), 0 2 (t), £(t), r|(t)}. That inertia coupling 
effect is manifested in the highly non-sparse structure of the non-constant, 
position-dependent "mass-matrix" M(y) in Eq.(3) which has the form (it turns 
out that mjj = mjj(6i, 0 2 ), in general; see Appendix A of this chapter for 
explicit my expressions) 



We will now proceed to show how one can develop TI and AD pointing 
controllers for the HST model in Figure 2, based on the "exact" equations of 
motion Eqs. (3), (4). 
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accommodated by well-known, simple modifications of the results presented 
here. 

Design of the TI mode of control begins by identifying the particular 
2 nd order equation of motion that governs the HST main-body rotation 4>(t) in 
Figure 2. That equation can be obtained by first multiplying Eq.(3) by M' l (y) 
to obtain (the existence of M -1 is established in Appendix B of this chapter) 

y + M _1 (y)Dy+ M-!(y)Ky + M“!(y) f(y,y) = M-!(y)bu (5) 

and then reading off the first of the five equations represented by Eq.(5). The 
result is the 2 n & order, non-linear equation of motion 

$ + g (0^02,01,0244^ ) = h(qi>q2) u ( 6 ) 

where the precise expressions for the functions g( • ), h( • ) are given in Eqs. 
(C.2) and (C.3) of Appendix C of this chapter. Eq.(6) governs the actual 
angular motions <|)(t) of the HST main body as shown in Figure 2, for an 
arbitrary control torque input u(t). Note that the one term g( *,•,*) in Eq.(6) 
embodies all of the disturbing torques induced on the HST main body by the 
flapping motions of the solar arrays in Figure 2. 

At this point, it is useful to introduce the notion of an ideal-model for 
the desired controlled (closed-loop) motions of the pointing angle <(>(t), in the 
TI mode of control. For this purpose, we assume the desired transient 
behavior of <|)(t) in the TI mode is represented by the solutions of the specified 
(given) "ideal-model differential equation" 

(Ideal Model Eq.) $ + (2£<a n )<fr + (<o n 2 )<|> = 0 (7) 

where the values of the parameters (£>0, co n >0) are assumed specified. 
Thus, the control designer’s task is to design the control function u=u(?) in 
Eq.(6) to make Eq.(6) "look like" Eq.(7). In the idealistic case, where all 
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arguments in g( • ), h( • ) can be accurately measured in real-time, it is clear 
that the ideal choice for the TI controller u( • ) in Eq.(6) would be 

u = h-i(0 l5 0*) [ - k 2 <j>- M> + g( 0 b 0 2 , 0 1? 0 2 i,4,fi) ] (8,a) 

where 

ki = co n 2 ; k 2 = 2 £ co n . (8,b) 

However, for this study it is assumed that only <|>(t) can be measured in real- 
time. Thus we will seek a physically realizable approximation to Eq.(8) as 
follows. First, the term h(0i, 0 2 ) in Eq.(6) will be approximated by a fixed 
constant h that is selected to approximate the range of values of h(0j, 0 2 ) for 
a representative range of (0i, 0 2 ) values. (The determination of h is 
presented in Appendix B of this chapter.) Next, the function g( in 
Eq.(8) will be viewed and treated as an unknown, time-varying ’’disturbance" 
term g(t) defined by 

g(t) = g( e 1 (t),0 2 (t),0 1 (t),0 2 (t),(j>(t),4(t),f|(t) ) (9) 

where (0i(t), 0 2 (t),(|)(t), £(t), q(t)} are arbitrary solutions of the exact 
equations of motion Eq.(3). A typical time-plot of Eq.(9) is shown in Figure 
4, where it can be seen that the unknown "disturbance" g(t) is a smooth, well- 
behaved function that slowly meanders back and forth in an oscillating 
manner. We will now show how one can estimate the time-function g(t) in 
real-time, from real-time measurements of <|)(t). Thus, following the principles 
of Disturbance-Accommodating Control Theory [3, p. 412], we first represent 
the uncertain time-behavior of g(t) by a quadratic poivnomial-spline (See 
Eq.(2).) 

g(t) = Ci + C 2 t + C 3 t2 (10) 
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where the weighting "constants" Q in Eq.(10) are allowed to jump in a once- 
in-a-while, stepwise-constant manner, as described below Eq.(2). Next we 
introduce a state-variable model for g(t) by defining 

Zl=g(t) ; Z 2 = g(t) ; z 3 = g(t) (11) 

and observing that the zj obey the following differential equations 

Z!=z 2 + CTi(t) ; z 2 =z 3 + o 2 (t) ; i 3 = cr 3 (t ) (12) 

where the symbols cjj(t) in Eq.(12) denote totally unknown, unmeasurable, 
sparsely populated time-sequences of randomly-arriving, random-intensity 
impulses (Dirac impulses). Those impulses represent the "cause" of the once- 
in-a-while jumping of the Q values in Eq.(10). 

The foregoing arguments and approximations enable us to rewrite the 
exact equation of motion Eq.(6) for <|>(t) in the approximate and simplified 
form 

$= hu - w(t) ; h«h ; w(t) =g(t) (13) 

where, according to Eq.(ll), the "disturbance" term w(t) in Eq.(13) is given 
by 

w(t) =Zi(t) (14) 

and where zj(t) is (approximately) governed by the set of impulse-forced 
differential equations Eq.(12). The TI control law (8) corresponding to (13) 
thus becomes the constant-gain linear control law 

u = h- 1 [-k 2 $-ki<j) + Zi(t) ] (15) 

where k lr k 2 are defined as in (8b) and where, following standard DAC 

A 

techniques [3, p. 430], <j>(t) and zj(t) denote on-line, real-time estimates of 
<j>(t) and zi(t) obtained from a composite plant/disturbance state-observer 
(Kalman filter), based on the approximate plant-model Eq.(13), and the 
"disturbance" model, Eq.(12). That composite state-observer operates on the 
plant output measurement <j)(t) and plant control input u(t) to produce accurate 
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real-time estimates of the state-variables {xj = <)>, X 2 = § } for the actual plant 
Eq.(6), as well as estimates of the state-variables { Z\ = g(t), z 2 = g(t), 
z 3 = £(*) } for the "disturbance" term w(t) in Eq.(13). The design of the 
composite state-observer for this TI mode of control is described in Appendix 
C of this chapter. 

Substituting the control expression Eqs. (15), (8b) into the actual 
equation of motion, Eq.(6), yields the actual, closed-loop equation of motion 


for <(>(t) in the TI mode as 

$+g(-,- > -) = h(-,-)h- 1 [-k 2 |-k,(|>+z 1 (t)] (16) 

which can be rearranged to read 

$ + h h" 1 [ ^©n)! + (cG n 2 ) 4> ] = -g( • , * ) + h h- 1 zi(t) (17) 

Comparison of Eq.(17) with Eq.(7) shows that the desired ideal-model 
behavior of (j)(t) in the TI mode of control will be realized, provided the 
assumed approximations 

h(- , ^fr 1 « 1 ; | « <j>(t) (18a) 

z i(t) « g( e 1 (t),e 2 (t),e 1 (t),0 2 (t),(i)(t)4(t),f|(t) ) (i8b) 


are achieved by the aforementioned design procedure. This completes the 
design of the HST controller for the TI mode of control. 

Design of the AD Controller 

The TI mode of control automatically isolates the HST main-body from 
any disturbance torques induced by the flapping solar arrays and is therefore 
proposed as the normal control mode for HST pointing experiments. On the 
other hand, one can imagine exceptional situations in which the prompt 
damping of excessive solar array oscillations becomes a primary concern. 
For those situations, the Array Damping (AD) mode of control provides a 
means for actively augmenting the natural damping of the solar array 
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oscillations by carefully controlled and orchestrated dynamic motions of the 
HST main body. To achieve this active augmentation to the solar arrays’ 
natural structural damping effects, the AD controller must maneuver the HST 
main-body motions <|)(t) back and forth, so as to induce controlled, time- 
varying, "bending torques" [at the body/solar array attachment point 
(structural interface)] that are strategically timed to dissipate the energy of the 
oscillating solar arrays themselves. It is rather surprising that this seemingly 
complex, intricate control task can be accomplished by a simple, constant- 
gain, linear controller, as we shall now show. 

The design of the AD mode of control begins by identifying the 
particular sub-set of the complete equations of motion, Eq. (3), that governs 
the motions <j>(t), 0i(t), © 2 (t). That particular sub-set consists of the first three 
2 nd order differential equations in Eq.(5), which can be separated out from 
(3) and be written as 

f + D[y)y + K[y)y + f(y,y) = b[y)u (19) 

where y = ( <j>, 0j, 02) and where D(y), K(y), f(y,y), b(y) can be determined 
from M^D, M -1 K, etc. in Eq. (5). Next, one must linearize Eq. (19) in the 
neighborhood of the anticipated "operating-point" values of {<j)(t),0i(t),0 2 (t)}. 
For our problem those operating point values are 

4> =<}>= 0 ; $1 = ^=0 ; 0 2 = 0 2 =0 (20) 
Linearizing Eq. (19) in the neighborhood of Eq. (20), using standard 
procedures, then yields Eq. (19) in the linear, constant-coefficient form [Note: 

f(y,y) = 0 at <j) -<j>= 0j = = 0 2 = 0 2 ~ 0] 

y + Dy + Ky = bu (21) 
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where { D,K,b } denote the values of {D(y),K(y),b(y)} at the operating 
values Eq.(20). In the AD control mode the objective is to simultaneously 
achieve the three stabilization conditions 


<(>(t) — > 0 ; 0 i(t) — ► 0 ; 02 (t) — > 0 ( 22 ) 

For this purpose we first seek u(-) in the (idealistic) constant-gain linear state 
feedback form 


u=kj 


(A 


- k cl (j) + k c2 0j + ^c 3^2 + ^c6^ : 


(23) 


where the six elements of the control gain-vector kc = (k^, k^) are 

chosen to place the corresponding six (closed-loop) eigenvalues 
of Eq. (21) at suitably "stable" locations in the complex plane. An effective 
procedure for computing the gains {k^, for any given set of desired 

eigenvalues {Xj, is outlined in Appendix C. For our simulation 

studies, we selected the six >4 to consist of three identical pairs of stable, 
complex-conjugate roots defined by 

; j=V-T , (24) 


where (£> 0 , © n > 0 ) are chosen to yield satisfactory closed-loop response 
characteristics for the AD mode. 

The final step in the design of the AD control mode is to develop a 
state-observer (Kalman filter) that generates accurate real-time estimates of 

4 >(t), 0 j(t), 0 (t), 0 2 (t), 0 2 (t) from the available measurement <j)(t), as needed 


to implement the AD control law Eq. (23). For this purpose we developed a 
conventional full-order observer based on the linearized model Eq. (21). The 
details of that observer design are presented in Appendix C of this chapter. 

^ A ^ A C 1 

The outputs (j>,0i,0i f 0 2 ,0 2 of that observer are used in place of their 
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idealized counterparts in Eq. (23) for implementation of the AD control mode. 
This completes the design of the HST controller for the AD mode of control. 

7. SUMMARY OF THE PROPOSED NEW DUAL-MODE 
DISTURBANCE-ACCOMMODATING POINTING CONTROLLER FOR 
THE HST 

The new HST pointing controller proposed in this study, for the single- 
axis model in Figure 2, consists of two distinct controller-algorithm 
configurations, one for the total isolation (TI) mode of control and one for the 
array damping (AD) mode of control. The TI controller is given by Eq.(15) 
and its associated state-observer [Appendix C, Eq. (C.13) ], where the 
parameters (£J , to n ) are assumed chosen to yield the desired quality of closed- 
loop regulation response <Kt) -* 0. The AD controller is given by Eq.(23) 
and its associated state observer [Appendix C, Eq. (C.26) ], where the design 
parameters ( £ ,o» n ) can be chosen to yield the desired "settling-time" for the 
actively augmented damping of the solar array oscillations. 

In both the TI and AD modes of control, the presumed real-time plant 
output-measurements, for the planar-motion configuration of Figure 2, consist 
of the HST main-body pointing angle <K0 only. If it should turn out that one 
can also accurately measure the rate <j>(t) in real-time, that additional 
measurement can be incorporated into the composite plant-state/disturbance- 
state observers associated with Eqs.(15) and (23); see Appendix C of this 
chapter. 

In any consideration of a dual-mode controller, the question of how to 
gracefully "switch" from one mode to another naturally arises. In principle, 
such mode switchings are usually accomplished by a slow ’’fading" from one 
mode to another, similar to the way one "fades" signals from right to left 
speakers in a stereo sound system. However, this fading procedure must be 
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done in such a way that one control mode does not tend to "fight” the other 
control mode during the fading/mixing process. In this project it was found 
that a "fade" to the zero control level, in one control mode, was necessary 
before attempting to fade-into the other control mode. Otherwise one ran the 
risk of "fighting" developing between control modes during fading, as alluded 
to above. A further investigation of this issue seems warranted. 

8. GENERIC PARAMETER VALUES FOR SIMULATION EXERCISES 
OF THE CLOSED-LOOP PLANAR-MOTION HST MODEL (FIGURE 
2 ) 

The numerical parameter values provided for the mathematical models 
supplied by NASA for this project were not appropriate for the alternative 
"first-principles," planar-motion dynamic model Eq. (3) we derived for our 
controller design procedure. Consequently, for our closed-loop simulation 
studies of the HST model in Figure 2, we chose a set of numerical parameter 
values, (for the various masses, inertias, length, etc. indicated in the 
configruration model of Figure 2 and in the analytical model Eq.(3) ) which 
seem to be representative of the relative scale of values associated with the 
real-life HST. As a matter of fact, since our "exact" model Eq.(3) is derived 
in symbolic mass, inertia etc. parameter terms, and our controller design 
procedure is likewise expressed in terms of those same symbolic terms, it is a 
simple matter to re-evaluate our controller expressions and our simulation 
results for any given numerical values of HST/Solar Array parameter terms. 
Accordingly, the following two sets of numerical parameter values for Figure 
2 were chosen for our simulation studies. 

For both the symmetric and asymmetric solar array cases 

loi=2.5m ; l 0 2=2.5m ; l g = 1.5m 
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lj = 4.8 m 
Mo = 10,5000 kg 
M! = 180 kg 
M 2 = 180 kg 


1 2 = 4.8 m 
J 0 = 10,5000 kgm 2 
Jx = 180 kgm 2 
J 2 = 180 kgm 2 


( 25 ) 


For the case of dynamically symmetric solar arrays 
At the 0.6 Hz. flex-mode. 

k n = 5.5 ; k 12 = 0.001 

k 2 i = 5.5 ; k 22 = 0.001 


At the 0.11 Hz. flex-mode. 
ku = 0.18 
k 2 i = 0.18 


k 12 = 0.001 
k 22 = 0.001 


For the case of dynamically asymmetric solar arrays 
At the 0.6 Hz. flex-mode. 

k n = 6.5 ; k 12 = 0.01 

k 21 = 5.5 ; k 22 = 0.001 


At the 0.1 1 Hz. flex-mode. 
k n = 0.3 

k 2 i = 0.18 


k 12 = 0.001 
k 22 = 0.001 


(26a) 


(26b) 


(26c) 


(26d) 


The set of HST controller parameter values chosen for use with Eqs.(25),(26) 
are: 
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For the TI control mode: 

C = 0.9 ; (Dn-1.0 (27) 

For the AD control mode 

C = 0.7 ; co n =1.0 (28) 

9. RESULTS FROM CLOSED-LOOP SIMULATION EXERCISES OF 
THE HST MODEL, EQ.(3) 

The "exact” non-linear model Eq.(3) associated with Figure 2 was 
simulated and exercised in a series of runs on a digital computer, using the 
two cases of plant parameter values shown in Eqs.(25),(26). In that series of 
simulation exercises, the control term u in Eq. (3) was implemented, first as 
the TI controller Eq. (15) and associated state-observer, and then as the AD 
controller Eq. (23) and its associated state-observer, using the corresponding 
controller parameter values shown in Eqs.(27), (28). To demonstrate the 
ultimate performance capabilities of the TI and AD control modes, control 
actuator saturation effects were not simulated in these concept-demonstration 
exercises . The simulated excitation of the solar array oscillations was 
accomplished by applying random, once-in-a-while, short-duration, high- 
intensity external torque pulses (simulated torque "impulses") to one or the 
other of the solar-array arms shown in Figure 2. The resulting angular 
motions <j>(t), 0i(t), 02© of the HST main-body and flapping solar arrays, for 
both the TI control mode and the AD control mode, and for both the 
symmetric and the asymmetric solar array cases, are shown in the series of 
time-plots presented in Figs. 5-7. The time-plots in Figure 5 clearly show the 
TI controllers' ability to quickly adapt to, and counteract, the time-varying 
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disturbance torques induced by the vibrating solar arrays and thereby 


THIS SPACE IS INTENTIONALLY LEFT BLANK 
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0.00 10.00 20.00 30.00 40.00 50.00 ^ 60.00 

Fig. 5a: Closed-Loop Response of <|)(t). 



0.00 10.00 20.00 30.00 40.00 50.00 ^ 60.00 

Fig. 5b: Closed-Loop Response of 0-j (t). 



0.00 10.00 20.00 30.00 40.00 50.00 ^ 60.00 

Fig. 5c: Closed-Loop Response of 02(t). 


Fig. 5: Closed-Loop Responses for the TI Control Mode with Dynamically Asymmetric SA's 
[0.1 1 Hz. SA oscillation triggered at t = 3.0 secs.] 
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Fig. 6: Closed-Loop Responses for the AD Control Mode with Dynamically Asymmetric SA's 
[0.6 Hz. SA oscillation triggered at t = 3.0 secs.] 
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effectively isolate the HST main-body from those disturbances, while 
maintaining a high-quality of pointing set-point regulation <j)(t) « 0. 

The closed-loop performances of <|>(t), 0i(t), 02© , corresponding to 
the AD control mode are shown in Figs. 6-7. Those time-plots demonstrate 
the ability of the AD controller to strategically maneuver the HST main-body 
motion <j>(t) so as to actively augment the natural structural damping in the 
solar arrays and thereby hasten the rate at which 0j(t) -> 0, 0 2 (t) -> 0. It is 
interesting to note that in the case of perfectly symmetric solar arrays. Figure 
7, the AD controller cannot distinguish which of the two solar arrays is 
causing the disturbing torques (i.e. a disturbance unobservability condition 
arises in the composite plant/disturbance state-observer used in the AD 
controller; see Appendix D of this chapter) and thus the AD controller cannot 
decide to which solar array it should direct the damping augmentation effort. 
Nevertheless, in that singular case, the AD controller automatically proceeds 
to maneuver the HST main-body so that the two solar arrays are coerced into 
flapping in an equal-but-opposite (birdwing-like) manner, 0j(t) = - 02 ©, so 
that the corresponding net torque disturbance induced to the HST main body 
is thereby damped to zero, while 0j(t) -* 0 and 02 © -» 0 according to their 
natural damping characteristics. 

When the two solar arrays are simulated as asymmetric , the AD 
controller succeeds in individually augmenting the natural damping of both 
0 i(t) -» 0 and 0 2 (t) -> 0 as evidenced in Figure 6. It is remarked that the AD 
mode of control, as presently developed, exhibits a relatively small domain of 
(closed-loop) stability and tends to be rather sensitive to system parameter 
variations. These undesirable features, which appear to be a consequence of 
the non-linear behavior of the terms M _1 (y)b and M -1 (y) f(y,y) in the basic 
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plant model Eq.(5), rather than a characteristic of the AD controller itself, 
constitute an important area for further research. 

10. A COMPUTER-ANIMATION VISUALIZATION OF THE 
SIMULATED HST CLOSED-LOOP RESPONSE DYNAMICS 

The unique performance characteristics associated with the TI and AD 
control modes are dramatically evidenced by viewing a computer-generated 
animat ion of the controlled movements of the 3-body HST model (main-body 
and two attached "arms" representing the two solar arrays) as depicted in 
Figure 2. For this purpose, the normal simulation data { <|>(t), 0i(t), 02© } 
was input into a specially prepared computer graphics program which 
repeatedly "draws" frames of Figure 2 for a sequence of closely spaced 
discrete-time values t^ , k = 0,1,2,3, ... . When that set of frames is displayed 
on the computer monitor in rapid sequence, one can "see" the three 
interconnected component bodies in Figure 2, "moving" in planar-motion just 
as they would in reality. The intricate maneuvers of <j>(t) one can thereby 
"see" being performed by the HST main-body, in the TI and AD control 
modes, are rather impressive. 

11. SUMMARY i CONCLUSIONS , AND RECOMMENDATIONS FOR 
FURTHER WORK 

In this chapter we have developed, and demonstrated by simulation 
exercises, a new dual-mode HST control concept for accommodating time- 
varying, persistent, uncertain torque disturbances (such as those due to flex- 
body "flapping" of the solar arrays) that act on the HST main-body and 
degrade precision pointing of the HST. The new controller consists of two 
distinct control modes. In the Total Isolation (TI) control mode (envisioned 
as the primary mode of control), the HST control torques automatically 
synchronize with, and adapt to, the time-varying disturbance torques (in an 
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equal-but-opposite sense) to effectively counteract and cancel the effects of 
those disturbance torques in real-time. As a result, the HST main-body is 
dynamically "isolated" from the disturbance torques, thereby permitting high- 
quality, precision pointing of the HST in the face of such disturbances. The 
system closed-loop performance using the TI control mode exhibits a high- 
degree of robustness. 

In the Array Damping (AD) control mode (likely to be used only 
occasionally, when not in a pointing experiment) the HST’s control torques 
are used to strategically maneuver the angular rotations of the HST main- 
body in a rocking, back-and-forth manner so as to induce an active damping 
(damping augmentation) effect to the natural damping of the solar array flex- 
body oscillations and thereby hasten the damping-out of those flex-body 
oscillations. When the AD control mode is used, the closed-loop system has 
the undesired feature of a relatively small domain of stability and relatively 
high sensitivities to system parameter variations, apparently due to the non- 
linear terms in the plant model, Eq. (5). A redesign of the AD controller to 
accommodate those non-linear terms should mitigate this feature. 

The qualitative performance capabilities of the proposed new HST 
dual-mode controller concept have been determined by closed-loop computer 
simulation exercises. Those simulations are based on the exact, non-linear 
equations of motion for a simplified, planar-motion, 3-body configuration- 
model (Figure 2) of the HST and its solar arrays using a generic set of mass, 
inertia, and other parameter values that are considered to be scalewise 
representative of the actual HST and its solar array disturbances. As 
indicated earlier in this chapter, the HST numerical parameter-values made 
available to us for this study (by NASA, Marshall) were not relevant to the 
"exact", single-axis nonlinear model (3) we derived and used here. The 
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closed-loop performance capabilities exhibited in our simulation studies of the 
simplified (single-axis) HST model are considered to be representative of the 
performance that can be achieved for the actual HST, using a full 3-axis 
(angular axes) controller designed by the same methodology employed here. 
The recommendations for further work, listed below, are directed at efforts 
that will validate this claim. 

Our recommendations for further work on the new HST pointing 
controller concept developed here are as follows: 

1. Generalize the configuration-model (Figure 2) and the associated exact 
equations of motion, to include arbitrary, 3-axis, coupled angular motions 
for both the HST main body and the attached solar arrays. Include 
multiple modes of solar array oscillations. 

2. Derive the 3-axis TI and AD controller equations corresponding to the 
generalized 3-axis model of HST. Generalize the AD controller equations 
to accommodate the non-linear terms in the HST model. 

3. Revise the 3-axis TI and AD controller equations, as needed, for 
implementation in digital control format, [4]. 

4. Demonstrate the closed-loop performance of the 3-axis HST model, and 
associated 3-axis TI/AD controllers, by computer simulation exercises, 
using actual HST controller actuator saturation levels and realistic values 
for masses, inertia, and other HST / solar array model parameters. 

5. Determine the most effective procedure for gracefully "fading” or 
"switching" from the TI control mode to the AD control mode, and vice 
versa. 
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Appendix A THE EXACT EQUATIONS OF MOTION FOR THE PLANAR-MOTION MODEL OF FIG. 2 

The exact equations of motion for the planar-motion model of Fig. 2 were created by the AUTOLEV 
program and output as a FORTRAN simulation. The pertinent definitions in the source code were then 
copied to MATHEMATICA and reduced to generate the equations shown below. Computer simulations of 
these equations were performed with the original AUTOLEV generated code to preclude the introduction 
of copying errors. 

For<£(*) 

— [Jo+Ji 4- 4 (& + A 2 *f l a 2 4- 2/i(/ 0 i cos -f l g sin 6\))Mi 

4(^02 4 A 2 4 l g 2 4 2A(A2 cos Bn — l g sin $ 2 )) M 2 ] <j> 

— [•A 4 (A 4 AiA cos$i 4 l\l g — [«A> 4 (l 2 4 A 2 A COS Bn 4 hlg sin -B-^M-i] Bn 

+ [(A + A sin^i)Mi 4 (l g - A sin #2^2] f 4 [— (Ai 4 A cos0i)Mi 4 (/02 4 A cos #2^2] 77 ( A. 1 ) 
= — u 4 <^/oiMi 4 r)<f>l g M\ — <p^l 02 M 2 4 r)j>l g Mn 

+ (£>£ 4 2<j>6il g 4 Q\lg)hM\ cos B\ 4 ( r}<j> — Ai — 0 2 Ai)AMi sin B\ 

— (j>£ 4 2(f>B 2 l g 4 8 2 l g )hM 2 cos Bn — ( 77 ^ 4 2 <fiB 2 /o 2 4 B 2 lc 2 )l 2 M 2 sm Bn 

For 0 2 (f) 

— 4 cosBi 4 hl g M\ sin ^1 ] 0 — [ J\ 4 / 2 Mi] 6 >i 4 [ 0]#2 

4[AAA sintfijf - [A Mi cos^i]r 7 (A. 2 ) 

= — [Ji 4 l\M\][k u Bi 4 Aroi^i] 4 AMi[(0£ — <^ 2 / 5 )cos^i 4 (<j>r) 4 0 2 Ai)sin#i] 

For & 2 (t) 

-[J- 2 +l 2 M 2 4 A2AM2 cos Bn — l 2 l g M 2 sin B 2 }4> 4 [O]0i — [J 2 4 / 2 M2]6A 

— [l 2 Mn sin Bn]£ 4 [InMn cos B 2 ]ij (A. 3) 

= — [ J 2 4 l 2 M 2 \[k 2 \B 2 4 ^ 22 ^ 2 ] 4 l 2 M 2 [(— $£ 4 <j>~l g ) cos B 2 4 (— <prj 4 <j> 2 lo 2 ) sin B 2 ] 

For £(t) [translation of m 0 c.g. in ^-direction] 

[Mi(lg 4 A sin B \ ) 4 M 2 (l g - In sin 0 2 )] i 4 [A Mi sin^]^ - [AM 2 sin02]0 2 - [M 0 4 Mi -fM 2 ]£4 [Oji; 

~ ~ [ifo(M 0 4 M\ 4 M 2 ) 4 4 >~(Iq\M\ 4 A2 Mo ) 4 

. . ... (A.4) 

A Mi (4" 4 2(pB\ 4 B\) cos^i — liM\(4>~ 4 2<j>Bn 4 0 2 ) cos # 2 ] 
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For rj(t) [translation of m 0 c.g. in ^-direction] 

— Mi/i[cos^i — cos 02 ]4> — [l\M\ cos0i]0i + {l\M\ cos0 2 ]02 4- [0]£ — [Mo 4- M\ 4- 

= [^(M 0 4- Mi + M 2 ) - 4> 2 l g (M x + Mi) - (A.5) 

1 1 Mi (<f> 2 4- 2<$0i 4" ) sin 6\ -+■ l2M2(4*~ -f- 2<£02 4" @ 2 ) sin ^ 2 ] 

State Variables and the Equations of Motion. The equations of motion (a.i)-(a. 5) can be repre- 
sented by the standard second-order matric-vector equation for the planar-motion HST configuration in 
Fig. 2: 


My 4- Dy 4- Ky = bw 4- /(y, y) 


where 


!/=[<£, 01, 02, rj\ T 

and the 25 elements [M] fJ of m are given explicitly by 

= — [0o 4- J\ 4- 0? 4* {^oi 4" ^i 2 4- l g ~ 4- 2/i (/qi cos 6\ 4- lg sin 0i) } M\ 

4* {^02 4" 4* lg~ 4- 2 / 2 (/02 cos 02 — lg sin 02) } Mi] 

[M]i, 2 = ~ [0i 4* 4* ^ 01^1 cos 81 -1- l\l g sin 8 \)M\] 

Mi* = ~ [*^2 4* (^2 4* loih cos 82 — hlg sin 6 i)M 2 \ 

[Mji^ = [( l g 4- h sin )M\ 4- (l g — h sin 02 )Mi] 

[Mji^ = [~(/oi 4- /1 cos 9\ )M\ 4- (/02 4- h cos 8 i)M 2 ] 

[M] 2 ,i = — (0i 4- l\M\ 4- Iq\1\M\ cos + l\l g M\ sin0i) 

[M]2,2 = — (0i 4- l\M\ ) 

[M] 2( 3 =0 

[M] 2,4 =l\M\ sin 0i 

[M]2,5 = - hMi cos 0i 

[M] 3i i — — (02 4* I0M2 4- I02I2M2 cos 0o — lilgMi sin 02) 

[M] 3i2 =0 

[M] 3i 3 = — (02 4" I 2 M 2 ) 


(A.6) 
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[M] 3 i4 = — l 2 M 2 sin 0 2 

[M]3 5 =/ 2 .M 2 COS 82 

[M]4,i = ( l g + l\ sin 6 \ ) 4- M 2 (/ 5 — I 2 sin # 2 )] 

[M] 4i2 =l\M\ sin 0i 
[M] 4t3 = — / 2 M 2 sin 0 2 
[M] 4 _ 4 = — (Mo 4* Mi + M 2 ) 

[M] 4i5 =0 

[Mjs,! = — Mi /1 (cos 0i — cos0 2 ) 

[M] 5 )2 = — l\M\ cos 0i 
—l\Mi COS 0 2 

[M]5 f4 =0 

[^1)5,5 — — (Mo + Mi + M 2 ) 

and where 

D =Diag [0, (Ji 4- /iAfi)i 12 , (J 2 + / 2 M 2 )fc 2 2 5 0, 0] 

K=Diag[0 , (Jj + /?Mi)A*n, (J 2 + /|M 2 )£ 2 1, 0, 0] 
b = [l, 0, 0, 0, 0] T 

The elements /,• of the 5-vector f( y , y ) are given explicitly by 

[f(y,y)]i = ^foiMi 4 - — 4>ilo 2 M 2 4 - r}4>l g M 2 

+ 4- c l(j) 6 ilg + 0jl g )l 1 Mi cos 0i 4- ( t)4> ~ 2<£0i/oi — ^2 ^01) fi Mi sin 0i 

- (4>i 4- 200 2 ^ 4- 0 llg)l 2 M 2 cos 0 2 - (t)0 - 1- 200 2 / O2 4- 0 2 / 02 )l 2 M 2 sin 0 2 
[f{y, y)]2 =/i Ml [(0f - 4> 2 l g ) cos 01 + (077 4- 0 2 /oi ) sin 0 X ] 

[f(y> y)]3 =hM 2 [(-<j>t + <j) 2 lg) cos 0 2 4 - (-077 + <j> 2 l 02 ) sin 0 2 ] 

[f{y,y)U = - [r)i(Mo 4- Ml 4- M 2 ) + 0 2 (/oi A/i + l 02 M 2 ) + 

l\Mi{(fr 4 - 200i 4 - 0[)cos0i — /iMi(0 2 4- 2 00 2 4- 0 2 )cos0 2 ] 
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[f(y i y)]5 — + M\ + m 2 ) — 4> 2 i g (Mi + M 2 ) — 

hMi(<j> 2 + 2 j>6i + 0 2 )sin0i + / 2 M 2 (<£ 2 + 2<£0 2 + 0^)sin0 2 ] 

The term Diag denotes a diagonal matrix with main diagonal shown. It should also be noted that the m 
matrix is strictly a function of e x and 0 2 (the matrices d,k are constants for the model (a.i)-(a. 5 )). The 
ten state-variables for ( a.6) can be chosen as (*! , . . . , x 10 ) = ( y , y ). 

Appendix B EXISTENCE OF M" 1 AND DETERMINATION OF A 

The existence of the inverse of the mass matrix, m, and of the scalar h = [M- 1 ]! 1 over the range 
of 0 i, $2 values of interest is easily established if the numerical values of the inertia, mass, and lengths 
in the planar-motion model are set. Using the parameter values (25) - (28) the determinant of m, for a 
representative range of 0i , 0 2 values, varies as shown in Figures B-1 . The variations in the h value for 
those same parameter values and the same range of e u o 2 values are as shown in Figure B-2. It is clear 
from Figures B-1 and B-2 that det . m and h are well-defined and strictly negative for the indicated range of 
0i, 0 2 values. The chosen approximating constant values for (m, h) are denoted by (M, A). 

For the linearization of (19) and design of the controller and observer gains, the representative constant 
values for m, m -1 and A were determined to be: 


1000M -1 



-30854.4 

-6487.2 

-6487.2 540 

i 0 





-6487.2 

-4327.2 

0 

0 

-864 




M = 

-6487.2 

0 

-4327.2 

0 

864 



(B.l) 


540 

0 

0 -10860 0 





0 

-864 

864 

0 

-10860J 




' -0.087898618 0.1317748 

0.1317748 


-0.0043706495 

0 


0.1317748 -0.43244034 

-0.19376112 


0.0065523381 

0.018988844 


0.1317748 -0.19376112 

-0.43244034 


0.0065523381 

-0.018988844 


-0.0043706495 0.0065523381 

0.0065523381 


-0.092298356 

0 

(B.2) 

0 

0.018988844 

-0.018988844 


0 


-0.095102461 




A = -0.000087898618 





(B.3) 


These constant linearized values were arrived at by setting 0 ! and 0 2 to zero and calculating the resulting 
values of m, m- 1 and h. 
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Appendix C DESIGN OF COMPOSITE PLANT/DISTURBANCE STATE-OBSERVERS FOR THE TI 

AND AD CONTROL MODES 

With the assumption that m is non-singular (a valid assumption as shown in App. B) then ( a.6) can be 
rewritten as 


y — — - M -1 Ky + M -1 bu 4- M~ x f(y, y) (C.l) 

DESIGN OF THE TI CONTROLLER GAINS. The equation for <j> contained in (C.l) is 


i = [M 1 ] i tx - g 


(C.2) 


where 

g =[M 1 ]i, 2 («/i + l\M\ 4- ^ 12 ^ 1 ) 

+ [M- 1 ] lt3 ( J 2 + llM 2 ){k 2 l e 2 -f k 22 e 7 ) (C.3) 

-M~ l f(y,y) 

The symbolic, analytical expressions for [M“ 1 ]i f i l [M- 1 ]i T2 , [m- 1 ^ are extremely long and complicated 
and are not reproduced here. Linearizing this equation at e x = e 2 = o results in 


i~ku — g(t) (C.4) 

where the g term includes the linearization error caused by the inevitable variation in m: that is, 


^) = 5-([m- 1 ] 1>1 -[m- 1 ] 1i1 ) u (C.5) 

Next, we approximate the time-variation of §(t) % Zl by the m’th order polynomial spline (here m = 
l , 2, 3, ... is chosen by the designer). 


g(t) = ci 4- c 7 t 4- . . . 4- 


(C.6) 


which implies that the dynamics of g{t) is modeled (approximately) by the m’th order system [3;p.413] 


'V 

h 


0 Im-1 


*1 

h 

+ 

’ <7i(0 

(m) 

L*lJ 


0 0 


(m-l> 

L z \ J 

&m— l(t) 


where i m _i is the m - i order identity matrix and indicates the m’th derivative of Zl . The elements 
<n(t) denote sparse sequences of randomly arriving, random intensity Dirac delta functions (impulses), 
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accounting for the once-in-a-while stepwise changes in the coefficients a of the approximating polynomial 
spline (C.6). 

Utilizing the control separability property inherent in DAC design methodology [3;p.438] allows the HST 
control torque action, u, to be split into a set-point regulation part u, and a disturbance cancellation part 
u d as follows 

u — u s -f u d (C.8) 

where one defines u 3 as 


u s — h 1 {k\(p 4- k 2 Q) ; assuming <p 3p = 0 


(C.9) 


and defines u d as 


u d — h l z\ (C.10) 

Assuming g{t) == £r, A « A, etc. the resulting closed-loop equation for 4 is: 

<j> + k 2 d) 4- k\<j> — 0 (C.ll) 

where ki and k 2 can be chosen as in (8.6) to obtain the desired closed-loop dynamics for <f>. 

DESIGN OF THE Tl OBSERVER GAINS. The composite equation of motion for (C.2),(C.5) and (C.7) 
is 


(m) 

Z\ 


0 Im+l 
0 0 




' 0 • 


' 0 ' 

4> 

Zl 

+ 

[M]‘ i 
0 

u + 

0 

(m-1) 

. Zi . 


. 0 . 


Pm (0. 


1 T 

Setting ^ Zl ... ( "7 I) J , (C.12) can be written 

x — Ax 4- bu + <r 

A full-order state-observer for (C.i3) has the form 

x = Ax -f- bu -f k 0 C(y — x) 


(C.12) 


(C.13) 


(C.14) 


UAH / ECE Dept. / Dr. C. D. Johnson 


1A-9 


Chapter 1, Final Report 


Design of a New High-Performance Pointing Controller for the Hubble Space Telescope 

NAS8-38609 / D. O. 49 


where (•) denotes the observer’s estimate of (■) and the plant output measurement is 


y = Cx = 0 (C.15) 

The state-estimation error for the observer (C.i4) is 

e-x-x (C.16) 

which obeys the well-known dynamical equation (between impulses of a(t)) 


s = (A--k 0 C)£ (C.17) 

where k„ = [* x , * 2 , . . . , & m+2 ] T . 

Using (C.12), and assuming that the oniy output measurement is <f>, the observer error dynamics (C.17) 
reduces to 


€ = 


—k\ 


km - f 1 
l—k m+ 2 


‘m-fl 


(C.18) 


The characteristic polynomial of (C.i8) is computed to be 


s m+2 + k lS m+1 + . . . + k m+1 s + k m+2 = 0 (C.19) 

The roots of this polynomial may be set to any desired values via proper selection of the observer gains 

kt. 

If both <j> and 0 are output measurements, the structure of the full-order observer (C.14) should be 
modified by setting y = (0, j>) T and redefining c accordingly, see [ 3 ;p. 430 ]. 

DESIGN OF THE AD MODE CONTROLLER. We set D = m - 1 d, k = m _ 1 k, / = m y ), and 
define the state vector x as 


x=ti>, 01, 00, i, 01, 0i\ T (C.20) 

Then, owing to the form of the d and k matrices, the state vector equation of motion governing x(t) 
becomes 


X 


— Ax + bu -(- £b 


(C.21) 
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where 


9 = {6 l +6 2 )/\/ 2 (C.26) 

A COMPOSITE STATE-OBSERVER FOR THE AD MODE. A full-order state-observer for (C.22) has 
the form 


x = Ax 4- bu + k 0 (y - Cx) (C.27) 

where I denotes the estimates of x and the output measurements are denoted by 

y = Cx = 0 (C.28) 

The observer’s state estimation error is defined by 


€ — X — X 


which has the error dynamics: 


(C.29) 


£ = (A - kt,C)£; k 0 = [k u k 2 , . . . , ke) T (C.30) 

Using the Linearized Modei in Appendix C, and assuming that the only output measurement is <f>{t), 
(C.30) reduces to 


r-*i 

0 

0 

1 

0 

0 ■ 

—k 2 

0 

0 

0 

1 

0 

— &3 

0 

0 

0 

0 

1 

— &4 


- [K]l ,3 

0 

-@1,2 


—^5 

-[K] 2 ,2 

-[K ] 2t 3 

0 

— [D]2,2 

— [6)2,3 

-*e 

“[K] 3 .2 

— [K ] 3 ,3 

0 

— [6)3,2 

— [6)3,3. 


The eigenvalues of (C.31) can now be chosen to obtain the desired observer dynamics and the corre- 
sponding k 0 can be computed. 

If both $ and i are output measurements the full order observer (c.27) should be reconfigured by setting 
y = [<f>, i \ T , modifying c accordingly, and proceeding as before. 
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Singular Case . As shown in Appen. D, the equation for computing the observer gain k 0 becomes sin- 
gular when m is symmetric, k n = *21 and k 12 — &22 • In this case, disturbances caused by the Solar 
Arrays appear as a single disturbance with one dynamical characteristic. In that case the same observer 
design method outlined above is used but with the transformation of variable (c.24) and associated state 
reduction (C.26). 


Appendix D CONTROLLABILITY AND OBSERVABILITY ANALYSIS FOR THE CASE OF 

DYNAMICALLY SYMMETRIC SOLAR ARRAYS 


Observability Analysis. When the solar arrays in Fig. 2 are dynamically symmetric (i.e. k n = *21 and 
£12 = k 22 ). The linearized state equations (C. 22 ) for 4, e u and e 2 result in a full-order observer equation 
of the form: 


x =A 0 £ -f b 0 u + k 0 (y — £); 

x = [<f) e l $ 2 $ 6 2 } t ; (D. 1) 


y -Cx - <f> 


where 


and 


‘0 

0 

0 

1 

0 

0* 


'0' 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

1 


0 

0 

6 

b 

0 

c 

c 

: b 0 = 

1 

0 

d 

e 

0 

f 

9 


0 

0 

e 

d 

0 

9 

/. 


0 


(D.2) 


k 0 — k 2 icq k$ ^6 j 

Assuming that y = Cx = <i> and defining the observer estimation error for (D.i) to be 

e = x — x 


(D.3) 


(D.4) 


The dynamics of e{t) are governed by 


i - (A 0 - koC)£ 


(D.5) 
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where 


—k\ 

0 

0 

1 

0 

O' 

—k 2 

0 

0 

0 

1 

0 

—^3 

0 

0 

0 

0 

1 

—k 4 

6 

6 

0 

c 

c 

— &5 

d 

e 

0 

/ 

9 

~k 6 

e 

d 

0 

9 

f 


Equating det [5/ - (A 0 - k 0 c)] to a general sixth order polynomial in 5 yields 


det [s/ (Ao — k^C)] — s 6 + P5S 5 + p 4 s 4 -f- p 3 s 3 + P 2 $~ ■+■ pis -f - Po (D.7) 

which allows the Jfe.’s to be determined. The result is: 

Ps = - 2/ + ki 

p 4 =(-2d + f 2 -g 2 )-2fk l + k 4 

p 3 =2 (df - eg) + (/ 2 - 2d- g 2 )k 1 + b{k 2 + k 3 ) - 2 fk 4 + c(k 5 -I- k 6 ) 

P2 ={d 2 - e 2 ) + 2 (df - eg)ki 4- ( cd + ce - 2bf)(k 2 + k 3 ) + (/ 2 - 2d - g 2 )k 4 

(D.8) 

4* (b — cf 4- cg)(k 5 4- &6) 

Pl =(d 2 - e 2 )^! + [b{e -d+f 2 -g 2 ) + c(g -f)(d + e)j {k 2 + *3) 4- 2(d/ - eff)* 4 
4- [&(<7 — /) 4- c(e — cf)] (£5 + k 6 ) 

Po = [b{d - e)(f + g) + c{e 2 - d 2 )) {k 2 + k 3 ) 4- {dr - e 2 )k 4 + b(e - d){k 5 4- k 6 ) 

As shown above, the six equations for Pi have only four variables: k 1 , {k 2 4- k 3 ), k 4 , and (k 5 + k 6 ). This is 
due to a loss of complete observability of the pair (a 0 , c) when the solar array flex modes shown in Fig. 
2 are dynamically symmetric. The observer gains are thus not generally solvable for the symmetric case. 
The set of equations (d. 8) however does suggest that in the dynamically symmetric solar array case the 
original state x in (D.i) should be reduced in dimension to read x = [<?, {$ 1 + $ 2 ), <j>, + e 2 )\ which 

will restore complete observability to the reduced system. 


Controllability Analysis. The controllability matrix for (c.22) in the case of symmetric solar arrays is 


[bo Aobo Agbo Agbo Agbo Agbo 


'0 1 0 0 0 0 ' 

0 0 0 0 0 0 

1 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 


(D.9) 


which obviously has rank =2 indicating that not all of the states of (C.22) are controllable when the two soiar 
arrays are dynamically symmetric. 
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Chapter 2 


FUNDAMENTAL CONCEPTS AND LIMITATIONS IN PRECISION 
POINTING AND TRACKING PROBLEMS 

Chapter Summary 

In this Chapter, we first describe the generic pointing and tracking 
problems in a general dynamical system/state-space context. Then, we 
analyze the information-theoretic aspects of the various uncertain signals in 
those problems, and establish some fundamental performance limitations 
those uncertainties induce, using various results and principles of modem 
control theory. It is shown that the introduction of "waveform models" for 
uncertain signals, leading to an extended-state formulation of pointing and 
tracking problems, is the most effective rational means of coping with those 
fundamental limitations. 

1. INTRODUCTION 

The design of pointing and tracking control systems constitutes one of 
the most important and widespread application areas for the field of control 
engineering. In spite of a long history of successful industrial applications, 
dating back at least to the early 1940’s, the fundamental scientific/theoretical 
aspects of pointing and tracking control problems are still not well-understood 
owing to the uncertain nature of several features of those problems. In 
particular, the time-varying pointing/tracking commands associated with such 
problems are typically not known a priori , but are only revealed on-line, in a 
real-time manner. Moreover, the internal and external disturbances that 
degrade the performance of pointing and tracking control systems are 
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inherently uncertain in nature and typically cannot be directly measured in 
real-time. These uncertainty features induce subtle performance limitations to 
the "solutions” of pointing and tracking control problems. The scientific 
aspects of those limitations have received very little attention in the literature. 

In this Chapter we examine the information-theoretic nature of the 
uncertainties inherent in pointing and tracking control problems, and identify 
the theoretical performance limitations induced by those uncertainties. It is 
shown that a novel "waveform model" representation of pointing/tracking 
commands, and of internal/ external disturbance actions, offers the most 
effective rational means of coping with those inherent performance 
limitations. 

2. POINTING AND TRACKING PROBLEMS; FUNDAMENTAL 
CONCEPTS 

The scientific consideration of pointing and tracking control problems 
necessitates a precise description of the essential conceptual features of the 
class of pointing and tracking problems to be considered. For this purpose 
we propose the (abstract) definitions given below. In both definitions, the 
dynamical system S being controlled is understood to consist of some 
physical hardware device (i.e. electro-mechanical-optical, etc.) that has a 
known, reliable mathematical model, a finite-dimensional state vector 
*=(*,, •••,*„) , r distinct control inputs w = ( u i9 u 2 ,— y u r ) (i.e. control 

motors, torquers, linear actuators, etc.) and p distinct disturbance inputs 
{wj(f),w 2 (/),.. ., >?,(/)} (i.e. wind gusts, target motions, coulomb frictions, 

structural distortions, flex-body vibrations, etc.). The outputs {y x ,y 2i — ,y M } of 
S consists of those features of S (positions, rates, accelerations, etc.) that 
can be measured by available sensors, in real-time. 
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Definition of the Pointing Control Problem - For purposes of this Chapter, 
the pointing control problem for a dynamical system S is understood to 
consist of the analytical design of the vector «(•) of control inputs 
i / = ( u n u 2 ,--,u r ) such that each of certain designated outputs y t (f) of S is 

forced to promptly attain, and steadfastiy remain at , a specified constant 
value (given set-point value) y i>sp in the face of an "arbitrary" initial-state 

x(t D ) of S and any transient or persistent vector w of disturbance-actions 
= ( w i (t),w 2 (t),---,w p (t)) that are likely to act on S . Generally speaking, 

the "constant" set-point values y iAp may abruptly change to some new 
constant value from time-to-time in an unknown, stepwise-constant manner. 

There are, of course, special applications of pointing control systems 
where it is required to "point" y i it) at one and only one, specific set-point 
y iJP , and therefore sudden, unexpected jumps in y iJp do not occur. To 

maintain the generality of our comments and characterizations we will 
hereafter disregard such special, restricted cases. 
Definition of the Trackine Control Problem - The tracking control problem 
for S consists of the analytical design of the vector of control inputs 
u = ( u i ,u 2 ,---,u r ) such that each one of certain designated outputs y(t) of S 

is forced to promptly coincide-with and thereafter faithfully follow (track with 
high-fidelity) an associated time-varying command input y lc (/), in the face of 

an "arbitrary" initial-state x(t o ) of S and any transient or persistent vector 
w(t) of disturbance actions that are likely to act 

on S. 

It is clear from these definitions that the general pointing problem is a 
special case of the tracking problem in which the time-varying "tracking- 
commands" y itC (t) become constant (set-point) commands y isp . 
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3 . FUNDAMENTAL LIMITATIONS IN POINTING AND TRACKING 
PROBLEMS: INFORMATION-THEORETIC ISSUES 

The primary sources of uncertainty, in the pointing and tracking problems 
defined above, he in: (i) the ways the set-point/tracking co mman ds {y isp (t) , 
y Uc (t)} and the disturbance inputs {w y (/)} naturally evolve in time and (ii) 
the way in which a pointing/tracking control system becomes aware of, and 
responds to, set-point/tracking commands and to the presence and nature of a 
disturbance input Wj(t) . 

3.1 Time-Evolution of Pointinz/Trackinz Commands 

In regard to the time-evolution of the set-point/tracking commands, it is 
typically the case that the time-behavior of {y Usp (t), y Uc (t)} , are not known to 
the control designer a priori. Moreover, the values of y irSp (t) and y i c (t) at 

time t are not known to the controller until time t . In fact, the underlying 
exogenous physical processes that determine how y iap {t) and y ic (t) evolve 
in time are typically such that at time t the future behavior of , 

y i c { t) , r> t , is not only not known, it is not determinable ! This latter fact 

is rather subtle and warrants further elaboration. Some typical time-plots of 
set-point and tracking commands, representing an after-the-fact recording of 
data, are shown in Figure 1 . The distinguishing feature of each plot, from our 
perspective, is that y isp {t) and y ic (t) have a smooth, well-behaved waveform 

pattern (trace) except at a few isolated random-like "transition times" t, 
where an abrupt change occurs in either the value of ( y itSp (t) , y t Jj)) and/or 
the value of some derivative of y itC {t). Between any two successive transition 
times t t , t M the plot of y itC {t) continues to be smooth and well-behaved, 

although not necessarily exhibiting the same pattern of behavior between each 
such pair of transition times. 
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If y iap (ti+e) is known, for a sufficiently small s > 0 then clearly the 
behavior of a set-point command y Usp {t) is accurately predictable between 
any two successive transition-times /, ,/ 1+1 . However, this seemingly 

deterministic character of y, i sp (t) is deceiving because the transition times 
t, themselves are determined by complex exogenous factors that are 

(typically) not quantifiable, a priori , and thus are inherently unknown, 
unpredictable and, in fact, are indeterminable, a priori. Therefore it is 
(typically) impossible to know or predict y itSp (t) beyond the (unknown) 

"next" transition time t M . 

A similar remark applies to tracking commands y ic (t) . Between any 
two successive transition times t t , t i+l the behavior of y i c (t) is 
conceptually predictable, at least over short sub-intervals of (/, ,/ f+1 ) , 
from knowledge of an initial record of y Uc {r) over a small interval t, 
Kr^iti + e )^ =small >0, since y ic (t) audits’ derivatives d k y t c (t) l dt k , k = 1,2,— , 
are all continuous in the open interval t { < t < t i+l . For instance, in that 
case one can argue that at each t a Taylor-series expansion of y Uc {t) 
should yield a non-zero sub-interval of convergence within the interval t f < t 
<t i+ 1 . 

The foregoing observations lead us to postulate the following prototype 
characterization of the class of realistic pointing commands y Usp {t) and 
tracking commands y f Jt) most often encountered in industrial applications. 
Prototype Characterization of Twical Pointing and Tracking Commands - 
The time-histories of typical pointing and tracking commands, illustrated in 

Figure 1, are characterized by a finite number of sequential, open, (not 
necessarily even), unknown intervals of time t, < t< t i+l , / = 0,1,2,*-, 

within which the functions y i sp (t) , y ic (t) and their derivatives, vary in a 
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smooth, well-defined manner. At the (unknown) times t-{ , t i+l ) the 

values of the functions y itSp {t),y Uc {t) , and/or one or more of their time- 

derivatives, experience simple jumps whose values are unknown and 
unpredictable, a priori . 

There are, undoubtedly, some exceptional cases one can cite in which it 
can be argued that this prototype characterization is not appropriate. 
However, we feel confident that this characterization is appropriate for the 
vast majority of real-life industrial applications involving pointing and 
tracking commands. 

Comparison with Random Process Characterizations - Some control 
engineers and signal-processing specialists, tend to view any signal that is not 
completely determinable (and completely known) as a "random" signal that 
can only be characterized by its (long-term) mean-value, variance, and other 
(higher) moments, as used in random process theories. Thus it is instructive 
to compare our prototype characterization of typical pointing and tracking 
commands given above to the stochastic characterization of random 
processes. 

The conceptual "signals" referred to in random-process theories as 
"white-noise" or "colored white-noise" [1] are mathematical constructs 
involving fictitious signals whose "values" go from "+oo to — oo" on every 
(arbitrarily small) positive interval of time. In particular, the Bode-Shannon 
realization of mathematical white noise consists of a densely populated time- 
sequence of Dirac-impulses having completely random arrival times and 
intensities, where the time-interval A t between any two successive impulse 
arrivals is arbitrarily small . Clearly, such a "signal" cannot have any 
waveform characteristics and, in fact, cannot exist physically. Moreover, if 
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such a signal is imagined to be fed into (input into) a physically existing 
linear dynamical system (a linear "coloring filter") the output response y(t) 

of that system (called "colored noise") will be devoid of any distinguishable 
waveform characteristics, since then some higher derivative of y(t) will 

always involve a "white noise" component. For instance, if one imagines the 

"coloring filter" to consist simply of a pure integrator so that 

y(t ) = £ (white noise)dz 

the response y(t) [called a Brownian or Wiener process] has no discernible 
waveform characteristics because its (formal) "derivative" dyldt is a densely 
populated sequence of randomly-arriving, random-intensity impulses (i.e. 
white-noise). 

It follows from these remarks that if, in the design of a tracking control 
system, the uncertain tracking commands y i c (t) are viewed (mathematically 

modeled) as "white-noise", or some colored hue of white-noise, the 

possibility of the design process accommodating any inherent waveform 
patterns that may characterize the actual physical commands y Uc (t) is lost 

forever. In particular, if one only models or specifies the long-term average 
statistical mean and variance (or power-spectral density) of y itC (t) in 

designing a tracking control system, one is (over)designing for an imagined 
class of commands (y l>c (/) } that includes outrageously rough, wildly-varying, 
continuously un-differentiable functions y i c (t) that would never occur in the 

actual physical application. Moreover, the long-term average statistical 
characterizations of mean, variance, power-spectral density, of y Uc {t) etc. 

may be completely irrelevant to the required short-term tracking performance 
of, say, a missile seeker tracking a maneuvering target. Such short-term 
tracking performance requirements demand that the seeker control system 
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react to the instantaneous waveform-pattern behavior actually exhibited by 
the target motions at each moment of time t — and such real-time waveform 
behavior information is never reflected in the long-term averages of mean, 
variance, etc. of target motions. The latter consideration is a manifestation of 
the old statistical joke about a person drowning while crossing a wide, placid 
stream whose average depth is 2 cm. 

In summary, high-performance tracking control systems operating in 
disturbance environments require real-time information about the 

instantaneous waveform characteristics of the particular tracking command 
y i<c {t) , and of the particular disturbance input w y (/) acting at each moment of 

time t . Statistical characterizations of long-term average mean, variance, 
power-spectral density, etc. do not (and cannot) provide such real-time 
waveform behavior information about tracking commands and disturbances. 
3.2 Characterization of Disturbances - The various disturbances «\(0 

that can act on the dynamical systems S involved in pointing and tracking 
control problems are almost always uncertain in nature and not accessible for 
direct measurement in real-time. Those disturbances arise from a variety of 
physical sources, some of which produce very erratic and capricious 
disturbance behavior. For example, fluid turbulence, radio-static, and similar 
"noisy" sources. Such disturbances have little, if any, discernible waveform 
characteristics and therefore, as a practical matter, can only be represented in 
terms of their averaged statistical properties. Thus, pointing and tracking 
control systems designed to cope with such noisy uncertain disturbances can 
only do so in a long-term averaged sense, such as minimizing the long-term 
average of the mean and/or the variance of the disturbance's effects on S . 
The design of such "good-on-the-average" controllers is addressed in the 
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subject of Stochastic Optimal Control Theory [2] and will not be further 
considered here. 

Fortunately, the vast majority of uncertain disturbances w.(t) that arise in 

realistic pointing and tracking problems are associated with physical sources 
that are not "noisy" in the sense of the preceding paragraph. Some co mm on 
examples of such uncertain disturbances are: coulomb friction and "stiction" 
effects, vibrations due to equipment motions, structural flexibility effects, out- 
of-balance and off-center effects, gyroscopic precession effects, gravity- 
gradient effects, crew motions, solar pressure effects, etc. In missile guidance 
problems the target's evasive maneuvers (jinking) appear as external 
"disturbances" in the guidance "error" equations. The time-history plots of 
such disturbances typically appear as shown in Figure 2, where it can be seen 
that, even though is uncertain and unpredictable, the time behavior of 

w.(t) tends to be smooth and well-behaved except at a few isolated, 
randomly occurring moments of time /, where the values of w ; (/) and/or 
some of the higher time-derivatives d k w j / dt k , k = 1 , 2 , 3 ,--, experience simple 

random-like jumps. In other words, the inherent uncertainty of such 
disturbances does not prevent them from being characterized in terms of 
certain trends or patterns of waveform behavior they naturally tend to exhibit- 
— between successive transition times (/ f ,/ i+1 ) . This feature of disturbances 
wfi) is recognized as being generically the same as that previously 

established for typical pointing and tracking commands; see Figure 1 . Thus 
we are led to the following characterization. 

Prototype Characterization of Typical Non-Noisv Disturbances in Pointing 
and Tracking Problems - The typical non-noisy disturbances {%(o} 

encountered in pointing and tracking problems are not known a priori , not 
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predictable, and not directly measurable. However the time behavior of each 
Wj(t) is characterized by a sequence of unknown, finite, open time-intervals 
during which w } .(t) and its derivatives are smooth, well-behaved functions of 
time. At the (unknown) ends of those time-intervals the value of w /(0 > 
and/or certain of the higher time-derivatives of w y (/), experience 
unpredictable, simple, random-like jumps. 

This completes our discussion of how typical set-point commands 
y isp (t) 9 tracking commands y Uc (t) and disturbances w f (t) naturally evolve 

in time. In the next section we will examine the performance consequences 
of the particular way in which pointing/tracking control systems become 
aware of, and respond to, pointing/tracking commands and the presence of 
disturbances. 

3.3 Poincelle’s Principle in Pointing and Tracking Control Problems - 
From the information-theoretic viewpoint, one of the primary contributors to 
performance limitations in pointing and tracking control problems is the fact 
that real-time changes which occur in the values of set-points or tracking 

commands at time t are not revealed to the control system until time t . 
Moreover, if some tracking command derivative d k y ic (t)/dt k , k = 1,2,—, 

abruptly changes at time t , the tracking control system typically will not 
become fully aware of, and properly react to, that derivative change until 

some elapsed-time after t , (here we are disregarding the unusual case in 
which y i e (t) and all of its time-derivatives can be accurately and directly 
measured in real-time). For example, if the value of dy i c (t)!dt abruptly 

changes at time t , a "type-2” (double-integral) tracking control system will 
not become fully aware of, and properly react to, that derivative change until 
the resulting growth in the tracking error 
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^(0=^, c (0-^(0 (1) 

has persisted for some (small) time-interval. Consequently, a momentary 
build-up of tracking-error (1) is necessary to initiate the proper controller 
response in such a situation. A typical time-plot illustrating this behavior is 
shown in Figure 3. 

The general fact that in a feedback control system some momentary build- 
up of error (1) is necessary to trigger the proper control response is known 
as Poincelle’s Principle. This principle serves to establish a limit on the 
quality of tracking fidelity that can be achieved in tracking problems with 
realistic, uncertain tracking commands of the type we are considering here. 

Poincelle’s Principle also applies to the effect of uncertain, unmeasurable 
disturbances wj(t) associated with pointing and tracking problems. In 

particular, the only way a pointing/tracking control system can become aware 
of, and properly react to, the presence of such disturbances is by virtue of the 
system output deviations (response errors) caused by those disturbances. 
Thus, some momentary build-up of pointing/tracking error (1) is necessary to 
trigger the proper control response to an unknown, unmeasurable disturbance. 

The only way one can overcome the pointing/tracking performance 
limitations imposed by Poincelle’s Principle is for the controller to have 

complete a priori knowledge of precisely how the "future’ 1 pointing and 
tracking commands {y Uip {t\y ic {t)\ and disturbances {w y (*)} vary 

with time, over the entire interval t a <t<T of the pointing/tracking control 
problem. We have already indicated that it is not realistic to expect such a 
priori information to be available in industrial pointing and tracking 
problems. However, it is interesting and instructive to consider how one 
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would design the "optimum” pointing/tracking control system for such an 
idealistic case. This subject is addressed in the next section. 

4.0 KALMAN'S THEORY FOR ABSOLUTE OPTIMAL CONTROL IN 
POINTING AND TRACKING PROBLEMS 

The performance limitations imposed by Poincelle’s Principle are a fact- 
of-life in virtually all realistic pointing/tracking problems. Nevertheless, it is 
interesting to investigate (theoretically) the nature of the optimum control 
system, and the degree of enhanced performance that would accrue, if in the 
formulation of a pointing/tracking control problem one could hypothetically 
have complete a priori information on precisely how the pointing/tracking 
commands and disturbances will vary with time over the entire fixture interval 
of control. Such a hypothetical problem was posed and solved by Kalman in 
a highly original 1963 paper [3]. In this section, we will summarize the 
essential features of Kalman's problem formulation and solution. 

4.1 Kalman's Pointinz/Trackinz Problem - The general pointing/tracking 
problem (called "servomechanism" problem) with disturbances, considered by 
Kalman in [3], was formulated for a very general class of multi-input/multi- 
output, possibly time-varying, linear dynamical systems S having the generic 
state-model [the symbols used below in (2) have been altered from those used 
in [3] to conform to current-day usage] 


[x = A(t)x+B(t)u + F(t)w(t)\ x(t 0 ) = x 0 
\y = C(t)x 


(2a,b) 


where * = (*,, •••,*„), */ = (*/,,— ,t/ r ), w = (w l9 - »,w p )andy = (y 1 ,-*, y m ) . The overall 
pointing/tracking performance in Kalman's problem was measured by the (still 
popular) "error-quadratic/control-quadratic" performance criterion J defined 
by the functional 
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J[u] = ^(r)&fr)+J[e r (0O('M0+« r (0^(0«(0^ (3a) 

e = (e,,-,e m ); e, -defined by (1), (3b) 

where the error/control weighting matrices are assumed symmetric 

and positive-definite for all t , and [ t a ,T ] denotes the specified time- 

interval over which the pointing/tracking control problem is defined. The 
problem is to find the control input w = (w, , -,w r ) to the dynamic system (2) 
that minimizes (3) for an a priori known pointing/tracking command input 

J'c(0=0i.(0,-,^(0). t.zt&T (4) 

and an a priori known disturbance input 

= h <t<T . (5) 

4.2 Kalman 's Solution - Kalman's solution to the optimum 

pointing/tracking control problem (2)-(5), with perfect knowledge of the 

"future," consists of the linear state-feedback control-law 

a* (/) = -Rr\t)B T {t)[P(t)x-h(t)] (6a) 

where the time-varying, symmetric positive-definite gain-matrix Pit) and the 
vector h{t) are determined, and stored for future use, by reverse-time solution 
[i.e. integration from t-T to t = t 0 ] of the following matric Riccati 
differential equation 

Pit) = -PA-A T P + PBR-WP - C T QC ; PiT) = C T SC (6b) 

and the coupled vector-matric differential equation 

hit) = i-A + BR-'B T P) T h + PFwit) - C T Qy c it) ; hiT) = C T Sy c iT) . (6c) 

From our perspective, the most significant feature of Kalman’s solution (6) is 
the fact that to carry-out the reverse-time solution of (6c), from t-T to 
t-t Q , requires complete a priori knowledge of y c it) and wit), over 
t a <t<T . In other words, after integrating (6) in reverse-time, at each t in 
forward-time playback (real-time use), the function hit) in the optimal control 
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law (6a) has embodied in it the optimal anticipatory control response that 
takes into account the (known) entire future behavior of y c (r) and w(r) , 

t<z<T . Thus, for example, if the tracking command or disturbance will 
have an abrupt jump in value at some ’’future time” t 0 <t x <T , the ’’all 

knowing” optimal control u°{t) in (6) "knows” that fact at the time t - 1 0 

and will begin optimally reacting to that future jump even before the jump 
occurs! Accordingly, we will hereafter refer to (6) as the absolute optimal 
control for the pointing/tracking control problem (2)-(5). 

The value of J~J° in (3) associated with a given initial state *&) of 

S in (2), and corresponding to use of the absolute optimal control (6) as the 
input to the dynamic system (2) , is the absolute minimal value of J 
achievable by any pointing/tracking control system, assuming the same 
x(t 0 ),y c (t) and w(t) ,t 0 <>t<>T . Kalman showed, in [3], that J° can be 

expressed as the sum of three terms: 

J° =x T (t 0 )P(t 0 )x(t 0 )-2xT{t 0 )h (t 0 )+u(t 0 ) (7a) 

where the scalar i >(/) is computed by solving, in reverse-time , the auxiliary 
scalar differential equation 

m = -yl {t)Qy c (t)+h J (t)B(t)R-Hl)B T U)h(t)+2h T (t)Fw(t); (7b) 

with the "initial condition" 

U (T) = yj (T)Sy c (T) (7c) 

The availability of J° enables one to explicitly compare the performance 
measure (3) for any practical pointing/tracking control system with the 
idealistic optimum performance (7 a) hypothetically obtainable if one could 
have complete future knowledge of y c (t ) , w(t) , a priori . In this way the 
fundamental performance limitations on (3), due to the practical inability to 
know the fixture behavior of y c (t) and w(t ) , can be rigorously quantified and 
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precisely measured for any candidate, physically realizable pointing/tracking 
control system. Although Kalman's solution was developed for the specific 
quadratic-type performance measure (3), his solution methodology can be 
applied to other pointing/tracking performance measures J , with similar 
results and conclusions; see [4], [5]. 

4.3 Apparent Failure of Bellman's Principle of Optimality - One of the 

cornerstones of modem optimal control theory for dynamical systems S is 
the assertion that the optimal control u°(t) at time / can always be 

expressed as a function of t and the system’s state x at time t . In other 
words 

, t 0 <t<T. (8) 

for some unknown function F (•, •) which is to be found by the analytical 
procedures of optimal control theory; see [6], [7]. 

Bellman popularized this assertion by elevating it to the status of a 
"principle” [8], which he called The Principle of Optimality. That principle 
can be paraphrased as follows: The optimal control at the present time t 
depends not on how the system S got into its present state xft) . but only on 
t and the value of the present state x(t) itself. In other words, the present 
state x(0 embodies all of the "past" that is optimizationally relevant to 
determining the optimal control at the present time t . 

A comparison of (6a) with (8) would seem to confirm Bellman's Principle 
of Optimality for the optimal pointing/tracking problem (2)-(5). However, 
when one realizes that the function h{t) in (6a) cannot be determined from the 
"present" data ( x(t),t ), but rather must be determined from the future 
behavior of y c (r) and w(r), t<r<T , it appears that the Principle of 
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Optimality fails in the case of pointing/tracking problems [note that the term 
Pit) in (6a) does not induce this characteristic]. 

In the next section of this Chapter we will show how the introduction of 
waveform-models. and associated state-models for uncertain 
pointing/tracking commands y i c {t) and uncertain disturbances w^t) enables 

the restoration of applicability of the Principle of Optimality to 
pointing/tracking control problems — at the expense of a higher-dimension 
state-vector and some inevitable loss of performance compared to the 
hypothetical, idealistic absolute optimal performance (7a). 


5. WAVEFORM-MODELS , STATE-MODELS , STATE-OBSERVERS 
AND A FUNDAMENTAL PRINCIPLE FOR ACCOMMODATING 
UNCERTAIN COMMANDS AND DISTURBANCES 

The effective solution of pointing and tracking control problems has 
traditionally been hampered by the lack of mathematical models that can 
effectively represent the kind of (non-noisy) uncertain pointing/tracking 
commands and disturbances encountered in realistic applications, i.e. model 
the prototype characterizations of commands and disturbances given earlier in 
Sections 3.1 and 3.2. In this section we will describe a relatively new form of 
mathematical model that is effective in representing those prototype 
characterizations of uncertain commands and disturbances. 

5.1 Waveform Models - The essential mathematical feature of the 
command and disturbance prototype characterizations given in Sections 3.1, 

3.2 is that they are '’analytic” between any two consecutive transition-times 
(/, ,f I+! ). An effective representation model of such functions is [here we use 
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the signal y i c (t) to demonstrate the ideas involved; analogous results apply 
to Wj(t)\ 

y,.c (0 = QJa (0 + Caf a (0 + — + Q*, /*, (0 (9) 

where the are totally unknown "constants" that are allowed to 

jump in value (at the transition times /,) in an unknown, once-in-a-while 
manner as shown in Figure 4 and the are completely 

known, smooth, well-behaved, analytic time-functions called "basis 
functions." The representation (9) is a generalized "spline-model" for the 
uncertain signal y Ue (t) and can be compared to similar representations used 


across spatial domains in "finite-element" techniques [9]. In practice the 
basis functions f } (t) are chosen by the pointing/tracking designer to reflect 

the variety of waveform behavior patterns actually exhibited by typical 
commands y ic (t). Thus, if the tracking command y ic (t) typically looks like 

that shown in Figure 5 the corresponding spline-model (9) would be chosen 
as 

^(0 = Q+Q<r* ; Le. /,(/) = 1; /„(/) = r» (10) 

Likewise, if an uncertain disturbance w } (t) typically looks like that shown in 


Figure 6 its spline-model representation (9) would be chosen as 

Wj ( t ) = C ;1 e~ 2i sin(3f ) 4- C j2 e~ 2t cos(3i) + C /3 (1 1 a) 

fix (0 = e ~ 2 ‘ sin(3/) ; = e- 2 ‘ cos(3 /) ; /,,(/) = 1 (1 lb) 

In some cases the typical waveform behavior of y t c {t) (and/or W j(0 ) is 
simply an uncertain, irregular meandering back-and-forth motion as shown in 
Figure 7. For such cases, the appropriate representation (9) is a "polynomial- 
spline" of the form 


y,.c (0 = Q\ + t+c i3 t 2 +‘ 


'Qm, * 




( 12 ) 
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Figure 4 - Stepwise-Constant Behavior of Weighting 
Coefficients C ;; in Waveform Model (9) 
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’Figure 7 - Plot of Uncertain Meandering Tracking Command that is 

Appropriately Represented by the Polynomial Spline Model (12) 
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where the value of m i = 1 , 2 , 3 ,- •• is chosen large enough to reflect the degree of 
back-and-forth motion exhibited by y ic (t) ; typically m i < 4 is satisfactory 

for industrial applications. Note that the case m i = l in (12) is the appropriate 

model for representing unknown, stepwise-constant pointing co mm ands 
y itSp 0). In (9)-(12) the unknown jumping of the "constants" C iS at the 

(unknown) transition times t i serves to represent the unknown jumping of 
y ic (t) , Wj(t) 9 and/or some of their time-derivatives, at those times. 

Tracking co mm ands and disturbances that admit a representation of the 
type (9) are said to posses "waveform-structure," and the representation (9) 
itself is known as a "waveform model." These concepts were introduced into 
the design of high-performance "disturbance control" systems in a series of 
papers beginning in 1968 [10]-[15]. The "information" encoded into the 
waveform model (9) is clearly non-statistical and is precisely the kind of 
dynamic, real-time command/disturbance waveform information needed in a 
pointing/tracking control problem to achieve maximum, physically realizable 
performance. 

5,2 Dynamic Models - In order for a pointing/tracking controller to access 
the information encoded into (9) it is necessary to re-format that information 
by converting the waveform-model (9) to an associated "dynamic-model", 
which is simply a homogeneous differential equation that (9) satisfies (for 
strictly constant, non-zero C i2 ). For this purpose one should seek the lowest- 

order of such an equation and, for mathematical convenience, it is prudent to 
select the set of basis functions {/j (/)} in (9) from among those that satisfy 

some linear (possibly time-varying) differential equation. Actually, it tums- 
out that for industrial applications it is usually adequate to choose basis 


UAH / ECE Dept / Dr. C. D. Johnson 


2-23 


Chapter 2, Final Report 





Design of a New High-Performance Pointing Controller for the Hubble Space Telescope 

NAS&38609 / D. O. 49 


functions f tj (/) in (9) that satisfy a constant-coefficient linear differential 
e quation . In the latter case one obtains a "dynamic-model" for (9) in the form 

d u, y itC {t)fdt Ui ^p^-\Xt)ldt^ x ^m^)ldt +M.W = 0 (13) 

where the constant coefficients and the order " u i " of (13), are 

completely independent of the C i; in (9) and are Ptffpisely .. determined, by the 
chosen basis-functions in (9). For example, in the case of 

(10) the associated dynamic-model (13) is 

X+3y, c = 0 ; (* ) = d()fdt (14) 

In the case (11a) the dynamic-model (13) is 

w] + 4 w t + 1 3 w i = 0 (15) 

The polynomial-spline waveform-model (12) leads to the m* order 

dynamic-model 

d m ‘y itC (t)/dr=0 (16) 

It is remarked that in some cases of multi-variable coordinated tracking 
c omman ds y c ~{y, c ,'",y mc ') (or multi-variable interacting disturbances) it may 
turn-out that the dynamic-model analogous to (13) consists of a coupled set of 
differential equations of the type (14)-(16). Note that (14)-(16) do not allow 
for (account for) the permitted, random-like, unknown once-in-a-while 
jumping of the weighting "constants" C fJ in (10)-(12). This deficiency will be 

corrected in the next subsection by introducing (symbolically) sparse 
sequences of unknown Dirac impulses as forcing functions in the "state- 
models" corresponding to (14)-(16). 

5.3 State-Models - The final step in converting (13) for use by a 
pointing/tracking control system consists of re-formatting each of (13) into an 
equivalent set of first-order (state-variable type) differential equations. Thus, 
for instance, by defining the "tracking command state-variables" for each 
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independent y i c {t) as: 

Command State- Variables = (c ( l , • ■ • • , c iUj ) 

(<4-i) 


( 17 ) 


the most general state-variable equation for a set of possibly coupled dynamic 

equations of the type (13) can be written as [12; p. 640] 

y c - [G](c) ; c=col. (c (1> |- • •|c ( ” ) ) = co/.(c, , • • -,c„); i>=j|U (18a) 

(<5) = [£](c)+M0 ; = (18b) 

where (G,E) are, respectively, mxv and vxv matrices that can be chosen 
from a family of canonical, completely observable pairs; see [12; p. 643]. 
The #.(/) in (18b) denote sparse sequences of unknown Dirac impulses that 
account for the random, once-in-a-while jumps that may occur in the C ij in 


(9). In a like manner, one can define "disturbance state-variables" for each 
independent Wj(t)as: 

Disturbance State- Variables z (0 = (z u , — , ) 


^j= w < 






'*Pt 


(Pi- 1 ) 
= W : 


(19) 


to obtain the following general state-variable equation for a set of (possibly 
coupled or coordinated) multi-variable disturbances w = (w, , — , w p ) 

(w) = [H]{z) ; z = col.(z il) \--]z (p) ) = coI.(z l ,--,z p ); p=£p i (20a) 
(z) = [D](z) + a(t) ; a=co!.(a u -",a p ) (20b) 

where, as in (18), the pair of matrices (H,D) can be chosen from a family of 
canonical, completely-observable pairs as explained in [13; pp. 405, 417] and 
the in (20b) denote sparse sequences of unknown Dirac impulses that 
account for the random, once-in-a-while jumping of the C i; that may occur in 


the disturbance counterpart of expression (9) , corresponding to the waveform 
model of A more general version of (20) is given later in (32). 

The particular choices (17), (19) of command and disturbance state- 
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variables are the most natural and convenient choices. However, as is well- 
known in the "state-theory" for dynamical systems, there are many other 

choices for state-variables that are equally valid [16; p. 170 ] — and some that 
are not valid; see [17 ]-[20]. The basis functions {/ y }, corresponding to the 
waveform models of y ic (t ), w^t) appear in (18), (19) as the eigenmodes of the 

matrices E,D respectfully. For instance, in the case of m independent 
stepwise-constant pointing commands y isp (t), which have only the one basis- 
function /,(/) = 1 each , one obtains E = 0, G-I , v-m in (18). Note that it 
is not required that (E, D) be stable in order that (18), (20) model realistic 
commands and disturbances; i.e., see (12), (16). 

It should be emphasized that, with the exception of the unknown, sparse, 
(symbolic) impulse-sequences /*(/), o{t) that have been added, the state- 
models (18), (20) are precisely the same as (and contain precisely the same 
"information" as) the dynamic model (13) and the disturbance counterpart of 
(13) (not shown). The virtue of the alternative format (18), (20) lies in the 
explicit recognition of y ije (t), and their time-derivatives, as individual, 
independent state-variables {c t j , z if } and in the fact that the existing theory 
of "state-observers" [Kalman filters], which we will be using in the sequel to 
estimate those state-variables in real-time, is formulated for general state- 
variable models of the type (18), (20). The assumption that the unknown 
impulse sequences #(*)» o>(0 in (18), (20) are sparsely populated is critical 
to the existence of waveform structure for y ic {t) and wj(t), respectively. In 

particular, as the impulses in the sequences #(/), a f (/) become more densely 
populated (time-spacing between successive impulse arrivals approaches 
zero) the sequences jug(t), a t (0 begin to behave like the Bode-Shannon 
"white-noise" described in Section 3.1 with the result that y itC (t), w : (t) in (18), 
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(20) lose all their waveform structure and (in those cases where E,D in (18), 
(20) happen to be stable matrices) begin to behave like conventional 
"colored-noises" of random process theory. 

5.4 A Fundamental Principle for Uncertainty Accommodation in 
Pointing/Tracking Control Problems - The introduction in (17), (18) and in 
(19), (20) of the overall pointing/tracking command state c- (c^,c 2 ,— ,c 0 ) and 
overall disturbance state z = (z : ,z 2 ,---,z p ) is the key technical idea that allows 

practical realization of a pointing/tracking control system which achieves as 
near absolute optimal performance levels as is physically, rationally possible. 
This result, referred to as the Principle of Optimal Disturbance 
Accommodation in [4], can be re-stated for the class of pointing/tracking 
control problems considered here as follows. 

The Command/Disturbance Uncertainty Accommodation 
Principle for Pointing/Tracking Control Problems 
Suppose the commands {y I C (0} and disturbances {h> ; (0} have 

known waveform-structures and state models in the sense of 
(9), (18), (20). Then, for a broad class of pomtmg/tracking 
performance measures J , and given dynamic systems S , the 
current values of the "states" {*(/), c(/), z{t)} embody sufficient 

real-time "information" about S and the uncertain dynamic 
behavior of {y 1>e (0} and {w,.(0} to allow a rational, physically 

realizable, scientific "optimal" choice for the real-time control 
input u(t) — even though the future behavior of y i c (t), w } .(t) is 

not predictable at the time t . 
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The rationale for this principle is as follows. The state models for 
x(/), c(t), z(t) are entirely deterministic and known, except for the sparse 
sequences o;.(r), #.(/) of unknown impulses which cause corresponding, once- 
in-a-while, jumps in (c(t), z{t)). Since the sparse random-like arrival-times 
and random-like intensities characterizing those impulse sequences are totally 
unknown (have no statistical structure), there is no rational, scientific method 
for anticipating, and reacting a priori to, their jump effect on z(t)) in an 
optimization problem. Thus, the ’’optimum" physically realizable, rational 
control strategy for such a situation is to choose */(/), at each t, to optimally 
react to the "current" values of *(/), c(t) 9 z(t) as if there will be no further 
jumps in c(t), z(t) in the future. This policy is equivalent to ignoring the 
presence of the oj(0»MO hnpulses altogether, and can be viewed as a 
pseudo-restoration of Bellman's Principle of Optimality for pointing/tracking 
problems, in the extended (x,c,z ) -state space. Of course, the ignored future 
jumps that will invariably and repeatedly occur in c(/), z(t) do, in fact, affect 
the absolute optimal choice of u(t), at each t. Thus our "optimum" physically 
realizable, rational control strategy, based on current values of x(/), c(t), z(t) 
only, does not preclude the possibility of attaining "better-than-optimal" [but 
never better than absolute optimal!] performance by non-rational methods 
such as by fortuitous guessing about [gambler's luck, speculating on] the 
future impulse arrival times and their intensities. However such gambling in 
function spaces carries high risks with low probabilities of success and is 
definitely not recommended for precision pointing and tracking problems. 

The "command/disturbance uncertainty accommodation principle" forms 
the foundation for the practical design of a new class of high-performance 
pointing/tracking control systems which must cope with uncertain tracking 
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c omman ds and disturbances having the prototype characterizations defined in 
sections 3.1, 3.2 and illustrated in Figs. 1-3. That design procedure will be 
outlined in Section 6. First, however, we need to indicate just how the states 
x(t), c(t) and z(t) can be accurately estimated in real-time. 

5.5 State-Observers ( Kalman Filters) for Estimation of (x(t), c(t), z(t )) in 
Real-Time - According to the uncertainty accommodation principle for 
co mmand s and disturbances, as stated in the previous section, optimal 
pomtmg/tracking control laws u° - JT(* *) should be sought in the extended- 

state feedback/feedforward format [compare with (8)] 

= f (x(t),c(t),z(t),t). (21) 

Thus, in addition to designing the "optimal" function JX—) in (21) it is 
necessary to develop a practical means of generating accurate estimates 
x(t),c(t\z(t) of the three state vectors x(t),c(t),z(t), in real-time. Fortunately, 
an effective means for accomplishing this latter task is provided by the 
existing Kalman Filtering Theory in the form of what is generically called a 
"state-observer." In this section, we briefly summarize the methodology 
whereby such state-observers can be designed to produce accurate real-time 
estimates x(t),c(t),z(t). We will focus attention here on the simplified cases 
where the pointing/tracking commands y Usp {t),y ic {t), and outputs ,y m of S, 

can be measured with negligible "sensor-noise". Otherwise one should use 
state-observers based on the full Kalman Filtering Theory, [21]. 

Observation of the Pointing Command State c = fo.c,.— ,c..) - As indicated 

near the end of section 5.3, the state-variables c : corresponding to a given, 
stepwise-constant pointing command = OV*<>W"'>W) are m in number 

(u~m) and are defined simply as 

<\=y^ ; ^ =y 2jp ;•••; c m =y mjp (22) 
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Thus, in the special case of stepwise-constant pointing-commands, no 
co mmand -state observer is required because the direct real-time measurement 
of the m signals y Usp (t) constitutes a direct, real-time ’'observation" 

(measurement) of the associated command state c- (c^,— ,c m ). 

State-Observers for the Trackine Command State c = (c„c 2 , -,c 0 ) 
Assuming that y c {t) is not a stepwise-constant command, as just discussed, 
and that an appropriate command state-model (18) has been developed, an 
accurate real-time estimate c(t) of the command state c(t) can be generated 
[between arrivals of the -impulses in (18b)] from real-time measurements of 
the /w-vector y c (t) by the following (continuous-time) data-processing 
algorithm (full-order state-observer) 

d(t) = Ec-K o [y c (0-Gc] (23) 

where the vxm "observer gain-matrix" K 0 is chosen by the designer to assure 
prompt decay to zero of the state estimation error 

e c (t) = c(t)-d(t) (24) 

between successive arrivals of the unknown impulses //,.(/) in (18b). It is easy 
to show that e c (t) in (24) is governed by the homogeneous linear equation 

e c = [E+K 0 G]e c (25) 

Thus, if E,G are constant, K 0 should be designed as a constant matrix that 
places all eigenvalues X Uo of [£ + K a G] in (25) sufficiently deep in the left-half 

plane. Since the pair (E,G) is, by construction, "completely observable," the 
latter is a standard task (pole-assignment) in modem control theory, [16]. 
The case of non-constant E and/or G is treated in a similar manner; see [12; 
p. 641]. 

The full-order (i/ 1 -order) state-observer (23) can be replaced by any one 
of a variety of reduced-order state-observers — with the expense of a 
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somewhat more complicated mathematical structure, compared to (23). One 
such observer, having the reduced-order (u-m), is described in [13; p. 487]. 
Composite State-Observers for Simultaneous Estimation of the System- 
State jc = (x. _) and Disturbance-State z = (z..-->.z r ) - Unlike the 

command-state c(t), which can be estimated via (23) from direct real-time 
measurements of the tracking command y c (t), the disturbance-state z(t) can 
(typically) only be estimated from measurements of the outputs y : (/),**• ,y„(0 
of the system S being controlled, (i.e. one typically cannot directly measure 
w(t) in (20)). If the equations of motion for the state x(t) of S can be reliably 
"linearized" about some nominal operating point (or operating condition) of 
S, the states (x(t),z(t)) of both S and w(t) can be simultaneously estimated in 
real-time by a (frill-order) "composite state-observer" as follows. 

Suppose the state x(/) = (^(/),^(/),— .^(Ojof S is governed by the 


known, linearized state-variable equations 

c T J*(0 = A (t)x+B(t)u(t)+F(t)w(t) 
L \y(0 = C(t)x 


(26a,b) 


where {A(t), B(t) y F(t), C(t)} are known matrices, and assume that an 


appropriate disturbance state-model (20) has been developed. Then, by 
incorporating (20) into (26) one obtains the following "composite 


plant/disturbance state-model" (hereafter we omit the permissible argument 


(t) on A, B, F, C, H, D to simplify notation) 


ViJ [o D \[z)lo\ 1 [a(/)> 

(27a) 


(27b) 


By defining x = coi.(x\z), the model (27) can be re-written in the more 


compact notation 
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i = Ax+Bu+B (28a) 

y = Cx (28b) 

where the meanings of A,B,C,a are clear from comparison of (28) with (27). 
A full-order [(n+p)th- order] state-observer for generating an accurate 
composite-state estimate l = co/.(x|£), between arrivals of the unknown pit) 

impulses in (27a), is given by the continuous-time data-processing algorithm 
[13; p.431], [22] 

I = Ai + Bu(t) - K a [y - Cl] . (29) 

The composite-state estimation error s = x-l associated with (29) is 
governed by the linear homogeneous equation 

t-[J+££]g (30) 

and thus the observer gain-matrix K a in (29) should be designed so that 
s(0~*0 promptly, between arrivals of the a; -impulses in (27a). The latter 
task necessitates that the composite pair (A,C) in (28) be "completely 
observable," [a technical condition that is typically satisfied in industrial 
applications but, strictly speaking, is not automatically guaranteed, even if 
both (A,C) and (A 77) are known to be, individually, completely observable]. 
The design of K 0 is otherwise a standard task in modem control, [22]. 

The state-observer algorithm (29) can easily be generalized to include: (i) 
systems S whose vector of outputs y(t) in (20b) is a (linear) function of u(t) 
and/or w(t) of the form 

yit) = C(t)x + E(t)u(t) + Uw(t) (3 1 ) 

and/or, (ii) disturbance sources such that the 
dynamic behavior of w(t) depends linearly on the system state x(t) and control 
u(t) as follows [compare with (20)] 


UAH / ECE Dept / Dr. C. D. Johnson 


2-32 


Chapter 2, Final Report 




Design of a New High-Performance Pointing Controller for the Hubble Space Telescope 

NAS8-38609/D. 0.49 


w(t) = H(t)z+L(i)x (32a) 

£ = D(t)z + M(t)x+N(t)u(t) + o(t) (32b) 

The case (31) can arise when the sensors, used to measure the system 
outputs y^t), include "accelerometers." The case (32) can arise from "state 
dependent" disturbances. Generalizations such as (31), (32), including 
additional "noise" terms, are addressed in [13; pp. 423, 424], [22], [23] and 
[24]. 

In addition to (29), various forms of reduced-order, composite-state 
observers for (28) can also be used to generate real-time estimates x(t)J(t). 

The details of one such observer are presented in [22]. When the 
measurements of the pointing/trackmg commands (y Usp (t),y ic (t)) and/or the 

system outputs y f (/) are corrupted with (additive) "sensor noise" the 
corresponding state-observers (Kalman Filters) have exactly the same form as 
(23), (29) with the exception that the observer gain-matrices (K 0 ,£ 0 ) in (23), 
(29), are then "optimally" determined by the (assumed known) noise statistics 
and by certain auxiliary equations (Riccati differential equations) as defined 
in Kalman Filtering Theory, [21]. 

In s ummar y, when measurement noises are negligible, one can obtain 
accurate, physically realizable, real-time estimates {x(/),c(/),£(0} °f the three 
state-vectors {x(t),c(t),z(t)} as needed for implementation of optimal 
pointing/tracking controls (21) by using the command-state observer 
algorithm (23) and the composite-state observer algorithm (29), provided the 
state-equations for S are linearizable as in (26). If measurement noises 
associated with y i>c (t),yXt) are not negligible, the (K 0 ,K 0 ) matrices in (23), 

(29) should be determined by the formal procedures of Kalman Filtering 
Theory. 
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6.0 FORMULATION AND SOLUTION OF POINTING/TRACKING 
CONTROL PROBLEMS IN THE EXTENDED (x, c, z) STATE-SPACE 

According to the command/disturbance uncertainty accommodation 
principle, and related arguments presented in Section 5.5, the design of high- 
performance, precision pointing/tracking control systems should be based on 
the real-time information embodied in the real-time values x(t), c(t) f z(t) of the 
three state-vectors x, c, z . That is, the designer should seek the 
pointing/tracking control law (algorithm) u = w(?) in the form 

"(0 c(Q, *(0, *) (33) 

where observer-produced estimates x(t),c(t),z(t) are used to implement (33). 

There are essentially two methodologies for designing extended-state 
pointing/tracking control laws of the type (33): (i) the optimal control method 
analogous to that used by Kalman for the problem (2 )-(5), and (ii) a more 
fundamental and direct method involving linear algebra and basic linear 
s tabilizati on techniques. In this section we will outline the essential features 
of both of those methodologies; further details may be found in [25], [26], 
[27]. 

6.1 An Extended-State Reformulation of Kalman's Optimal 

Pointinz/Trackinz Problem - Kalman’s absolute optimal, physically un- 
realizable, pointing/tracking problem (2)-(5) can be reformulated into a 
physically realizable, slightly less-than-optimal, problem by introducing the x, 
c, z state vectors, and the associated state models (2), (18), (20) with the 
sparse, unknown impulse-sequences ignored, [since there is no 

rational, scientific way (only "gambling” ways) to account for such uncertain 
o(t),ju(t) ia this, or any other, optimization procedure]. For this purpose, note 
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that the tracking error-vector e in (1),(3) can be rewritten in terms of (x, c, z) 
as follows 


e = y c -y 


Thus, 


= Gc-Cx 
= l-C\G\0] 


(x\ 

c 

Z J 


e-Ci\ C=[-C|G|o]; J = co/.(x|c|z) 


(34) 


Moreover, the plant, command and disturbance models (2), (18), (20) can be 
combined into the one composite i model 


A 

0 

FH 



(") 

r oi 

0 

E 

0 

(!) + 

0 

+ 


0 

0 

D 


0 


{ <**)) 


y - [c| o\ o] (f) 


(35a) 

(35b) 


and, using (34), the error/control quadratic performance criterion J in (3) can 
be rewritten in terms of x as follows 


J{u] = i T (T) 51(7)+^ ^ T QZ+u T Ru\dt\ 


S=C T SC 
Q = U T QU 


(36) 


The control w() that minimizes (36) for the composite 
plant/command/disturbance model (35) [with //(/) = o, o{t) = 0] can now be 
derived, without needing to know the future behavior of y c (t) and w(t), using 
standard procedures of optimal control theory. This has been done in [25] 
with the result that the ”less-than-absolute optimal," physically realizable, 
pointing/tracking control law (33), based on the extended (x,c,z) state-space, 
is given by 

Ho = -R-\t)B T (t)[K x (t)x+K^t)c+K :a (t)z] (37) 
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where the three gain-matrices (K^K^Kj) are determined by reverse-time 
solution of the following, unilaterally-coupled matric differential equations 
[25; p. 349] 

K x (t) = -K X A - A T K X + KJBR-'B t K x - C T QC ; K X (T) = C T SC (38a) 

JC(0 = -K xc E-A T K xc +K x BR~ l B T K xc +C T QG ; K xc (T) = -C T SG (38b) 
^ (/) = -KJFH- A t K b + - K b D ; K B (T) = 0 . (38c) 

Since the sub-absolute optimal problem (35),(36) is the same as Kalman’s 
problem (2)-(5), with the exception that in solving (35),(36) the ait), n{t) 
impulses have been ignored, it is interesting to compare the structures of the 
two control algorithms. Comparing (6a,b) with (37),(38a) it is clear that one 
always obtains the equivalence 

K x (t) = P(t) in(6a,b). (39) 

Moreover, if in fact the sparse impulse-sequences ait), //(/) in (35) are 
actually zero, then it can be shown [28] that the following additional 
equivalence obtains 

£„«)€+ K B (t)z = -h(t) in (6a,c) (40) 

Finally, still assuming o{t) = 0, /u(t) = o, it tums-out [25] that the extended 
state-space counterpart J° of the "optimal performance" expression J° in 

(7a), corresponding to (36), (37), (38), is given by 

> = ^ r ajK x 0 o )x(O+^(/ o )K c (/ o )c(O+ Ix^QK^iQciO 

+2[ xT (i 0 )K^O+c T (OKJt c )] z(t„) + z T (t 0 )K,(t„)z(t 0 ) (41) 

where the additional matrices K a {t), K c {t), K z (t) are determined by reverse- 
time solution of other matric differential equations similar to (38); see [25; 
Eq.(10)]. Note that, for fixed values of x(t o ), c(tj, z(t o ) and T, the "optimal" 
performance J 0 in (41) will always be larger (less optimal) than J° in (7 a) 
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when the o(t), fi(t) impulses in (35) are not identically zero . In that case the 
difference L defined by 

£-J°-J° (42) 

is a direct measure of the pointing/tracking performance loss due to the 
inability of the physically realizable, sub-absolute optimal control u°(t) in (37) 

to "see" the forthcoming, future impulse arrivals in a(0, m ( 0 • The practical 
implementation of (37) will, of course, require use of observer-generated real- 
time estimates *(/), c(t\ z(t) as developed in Section 5.5. 

An experimental comparison of Kalman’s ’’absolute optimal” performance 
j° in (7a) and the physically realizable sub-absolute ’’optimal" performance 
j° achieved using the extended-state optimal control design method (35-(38), 
based on simulation studies of a particular example, is presented in [28]. For 
the several cases studied there, it tums-out that the physically unrealizable 
absolute optimal performance J° is about 10-20% "better" (i.e. lower) than 
the physically realizable performance J°. 

We will now describe an alternative and more fundamental method for 
designing pointing/tracking control systems in the extended (x,c,z,) state- 
space. This alternative method does not involve optimal control ideas, and 
the introduction of a (potentially) "contrived" performance criterion J, but 
rather is based on simple linear-algebraic concepts and linear-stabilization 
techniques. 

6.2 An Algebraic-Stabilization Method for Pointing/T racking Control 
Design in the Extended (xx.z) State-Space - The optimal control method 
outlined in Section 6.1 is an effective control design methodology because of 
the wide-spread popularity of Linear/Quadratic-type methods among today’s 
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younger control designers and the wide-availability of a variety of (nearly 
automatic) CAD programs that are specifically tailored for executing 
Linear/Quadratic design procedures. On the other hand, such optimal control 
methods tend to relegate all conventional control design considerations (e.g. 
"setting-time," % overshoot, etc.) to the background and focus attention on 
mathematically minimizing the performance functional J. This introduces an 
extra layer of abstractions and mathematical complexities that can obscure the 
basic, underlying scientific features of the fundamental control problem being 
addressed. In this Section we will outline an alternative pointing/tracking 
control design method, for the extended (x,c,z) state-space, that is purely 
linear-algebraic and linear-stabilization in nature and that, like classical 
control methods, focuses on the essential, underlying features of 
pointing/tracking problems. Further details of this algebraic-stabilization 
method may be found in [26], [27 ], [28]. 

The algebraic-stabilization method of pointing/tracking control design for 
(2), (18), (20) [in (x,c,z)-space] begins by observing that the tracking error- 
vector e(t)-y c (t)-y(t) in (34) can become identically zero only if, for any 
arbitrarily given cjmmand-state c = (c 1 ,---,c y ) , there exists a system-state 
vector x such that 

Cx-Gc . (43) 

Satisfaction of the linear algebraic equation (43) is guaranteed if, and only if 
rank\C | G] = rank\C\ (44) 

which implies that G in (18) must allow the decomposition (factorization) 

G = C6 , for some (perhaps non - unique) matrix 6 .(45) 
The fundamental requirement (44) is referred to as the "trackability 
condition," [26; Thm. 2]. We hereafter assume (44) is satisfied and, for 
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simplicity, that A,B,F,C,H,D are all constant; see [26] for the more general 
time-varying case. Next, we re-write e as follows 


e = Gc-Cx 
= C(Oc~x) 

e ~Ce ss ; - {9c -x) - the " servo - state" vector , (46) 

It is clear from (46) that the basic requirement for achieving precision 
pointing and tracking, e(t) « zero , is to control the auxiliary " servo-state” 
vector e„ so that e u (t) is promptly regulated to, and thereafter remains in (or 
very near) , some subspace of the null-space OF [C] of the given matrix C, 
while remaining acceptably bounded, in norm, at all times. Recall that e a = 0 
is always a member of OF[C] ; however, other values of e u * 0 may also lie in 
the linear space of vectors that comprises OF[C\ , depending on the rank of 
C. 

The state equation governing the motions of e n {t) is readily computed to 
be 


= Ae u +(0E - A0)c-Bu-Fw(t) + 0p(t) (47) 

Thus, the control «(•) in (47) should be designed to promptly regulate 
> some subspace of OF [C], between successive arrivals of the unknown 
impulses in o(t), while keeping ||e„(/)|| acceptably bounded for all t . 
For this purpose, we will seek «(-) in the idealized, linear, extended-state 
feedback/feedforward form [26; eq. 36] 

u^S.x+S^+Tz (48) 

where the three constant gain-matrices S } , S 2 , r are to be designed as 
explained below. 

Substitution of (48) into (47) yields the ’’closed-loop" equation of motion 
for e^t) as 
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e a = (A +BS,)e u + (0E-A 6 - BS,0- BS 2 )c- (BT+ FH)z (49) 

Thus, to achieve prompt regulation of e a (t) -+S NC , S NC csome subspace of 
N[C ] , one should: 

(i) Design S l so that all solutions e„(t) of - (A+BS x )e a are strongly 

asymptotically stable to some S NC c N [C] and remain norm- 
bounded at an acceptable level, for all t ; see [29] for the theory of 
"subspace stabilization." (50) 

(ii) Design^ so that the vector (6E- A 0- BS ] 0- BS 2 ) c always lies in 

S„ c c N [C ] , for any c. (51) 

(iii) Design F so that the vector (Br+FH)z always hes in S NC c N [C ] , 

for any z . (52) 

(iv) Re-check satisfaction of the boundedness condition in (?) and revise 

choice of S NC and/or designs of S l1 S 2 ,r as needed; see [27]. (53) 

These design requirements necessitate the trial-and-error choice of a suitable 
subspace S NC c 0F[C] on which to "land" e a (t), but otherwise involve only 
linear-algebraic considerations that can be achieved by straightforward, 
systematic procedures as explained in [26], [27], [29]. 

An overall block-diagram showing the control (48) implemented on 
the system (2) [using full-order command-state and composite 
system/disturbance state-observers (23), (29)] is shown in generic format in 
Figure 8. A specific 2nd-order example which illustrates this algebraic- 
stabilization design method is worked-out in detail in [26; p. 3 1]. 

7.0 SUMMARY AND CONCLUSIONS - 

In this Chapter we have given a general description of a broad class of 
multi-input/multi-output pointing and tracking control problems encountered 
in industrial applications. The dynamic characteristics of the uncertain 
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Figure 8 - Generic Block-Diagram of Dynamic system (2) with Tracking- 
Command Model (18), Disturbance Model (20), and (x, y, z)- 
Based Tracking control (48) with State Observers (23), (29). 
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commands and disturbances associated with those pointing/tracking problems 
have been analyzed from the information-theoretic viewpoint and some 
fundamental pointing/tracking performance limitations induced by those 
uncertainties have been identified. It has further been shown that the solution 
of a theoretical "absolute optimal" pointing/tracking control problem without 
such uncertainties, due to Kalman, can be used to quantify and evaluate 
numerically the extent of those uncertainty-related performance limitations, in 
any given application. 

A relatively new form of mathematical model ("waveform-model"), which 
can simultaneously represent trend behavior and uncertainty in realistic 
pointing/tracking commands and disturbances, has been used to obtain a 
c omman d state-vector c(t) and disturbance state-vector z(t) 9 and associated 
state-variable models, that can effectively represent the characteristic 
dynamic time-behavior of uncertain commands and disturbances in 
pointing/trackmg problems. A fundamental principle, called the 
"co mman d/disturbance uncertainty accommodation principle," has been used 
to assert that an optimal, physically realizable, rational control law 
(algorithm) for a pointing/tracking control problem must be based on the real- 
time values of the current system-state x(t) , the current command-state c(t) 
and the current disturbance-state z(t). Two different control design 
methodologies, a variation of Kalman's optimal control method and a purely 
linear-algebraic method, for deriving such pointing/tracking controllers have 
been outlined. 

The fundamental principles and performance limitation ideas presented 
here, together with the extended-state (x,c,z) problem formulations and 
associated control design methodologies in Section 6.0, should enable the 
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attainment of maximum, physically realizable performance levels in a broad 
class of precision pointing/tracking control problems with uncertain 
commands and disturbances. 

The ideas presented here can be readily adapted to the design of digital 
(discrete-time) pointing/tracking control systems by following the procedures 
given in [14]. It is remarked that the use of a digital control system in 
pointing and tracking problems introduces an additional degree of 

performance limitation owing to the imposition of a specified sequence of 
discrete times t j at which the (otherwise "open-loop") discrete control actions 

are allowed to be updated. An effective technique for mitigating that 
additional performance limitation is described in [14; p. 311], [30]. 
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